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A DISCRETE LOG GAS, DISCRETE TOEPLITZ 
DETERMINANTS WITH FISHER-HARTWIG 
SINGULARITIES, AND GAUSSIAN MULTIPLICATIVE 

CHAOS 

CHRISTIAN WEBB 


Abstract. We consider a log-gas on a discretization of the unit circle. 
We prove that if the gas is not too dense, or the number of particles 
in the gas is not too large compared to the scale of the discretization, 
the absolute value of the characteristic polynomial can be described in 
terms of a Gaussian multiplicative chaos measure. This is done by ana¬ 
lyzing discrete Toeplitz determinants with Fisher-Hartwig singularities. 
In particular, we prove that if the gas is not too dense, the classical 
Fisher-Hartwig conjecture holds for the discrete Toeplitz determinant 
as well. Our analysis suggests that if the gas is any denser than this, 
the formula needs to be modified. 


1. Introduction 

It is a basic fact in random matrix theory that for a large class of random 
matrix models, the distribution of the eigenvalues can be described in terms 
of a log gas - namely a Gibbs measure of a system of particles experiencing 
an external confining potential and interacting with each other through a 
logarithmic repulsion. Notable examples of such models are the Gaussian 
ensembles, circular ensembles, and the Ginibre Ensemble - for more details, 
see e.g [25, 12], 

When considering the global (or macroscopic) scale, asymptotic proper¬ 
ties of the distribution of eigenvalues are often described in terms of linear 
statistics: if (Ai,..., Xn) are the eigenvalues, one considers some 

nice enough function g. The leading order behavior of such linear statistics 
are known to be described by a law of large numbers - the asymptotic dis¬ 
tribution is described in terms of a deterministic equilibrium measure (see 
e.g. [17]). The fluctuations around this equilibrium measure are known in 
many cases to be Gaussian (see e.g. [9, 17, 27]). 

In fact, these fluctuations in the linear statistics imply that the fluctua¬ 
tions in the distribution of the eigenvalues, e.g. on the level of the logarithm 
of the characteristic polynomial, are asymptotically described by the Gauss¬ 
ian Free Field (perhaps restricted to an interval, the unit circle, or the unit 
disk and with a suitable convention for the zero mode) - see e.g. [16, 13, 27]. 
Such fluctuations have been seen even on the level of discrete log-gases - for 
these, the underlying space where the log-gas lives on is a discrete set. Such 
models arise in for example combinatorial questions - see e.g. [18, 5]. 
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The Gaussian Free Field is a rough object (for an introduction to it, see 
e.g. [28]). In particular, its ’’geometric” properties are far from trivial. It 
has recently been discovered that the geometry of the Gaussian Free Field, 
and more generally log-correlated Gaussian fields, is best described by ’’ex¬ 
ponentiating” them into random measures through the theory of Gaussian 
multiplicative chaos going back to Kahane [20] (for a review see [26] and for 
a concise proof of existence and uniqueness, see [4]). This exponentiating 
has proven to be intimately related to the multifractal properties of the field 
(see e.g. [26]) and the extrema of it (see [10, 22]). These random measures 
also play a crucial role in the mathematical study of the random geometry 
of two-dimensional quantum gravity ([11]) and construction of conformally 
invariant random planar curves ([29, 1]). 

The way these measures are constructed is by regularizing the given log- 
correlated field into a function, exponentiating and normalizing this, and 
then passing to a limit where the regularization is removed. A natural 
question is then do all (reasonable) regularizations produce a multiplicative 
chaos measure in this way (and is the law of the limit the same for every 
regularization). Mostly the regularizations which are known to produce a 
chaos measure are trivial in the sense that the regularization imposes a 
martingale structure or something else which makes proving convergence 
quite simple. To the author’s knowledge, the only case, where a Gaussian 
multiplicative chaos measure has been proven to emerge from a ” non-trivial” 
regularization of a log-correlated field, is the characteristic polynomial of the 
circular unitary ensemble (see [32] based on conjectures in [14] and results 
in [7, 8, 6]). One of the goals of this note is to offer an example from another 
class of models - the characteristic polynomial of a discrete log-gas which is 
in fact a discretization of the law of the eigenvalues of the CUE. 

The main tool for analyzing our model is the realization that the relevant 
quantities can be described as discrete Toeplitz determinants with Fisher- 
Hartwig singularities. In our setting, a discrete Toeplitz determinant is one 
whose entries are given in terms of a discrete Fourier transform of a symbol 
instead of a continuum one. In [3], the asymptotics of such discrete Toeplitz 
determinants are studied in the case where the symbol has no singularities. 
Comparing to the continuum setting, it is a natural question then to consider 
what happens when the symbol can have Fisher-Hartwig singularities. A 
further goal of this note is to describe (in terms of the relationship between 
the number of particles in the gas compared to the number of lattice points 
in the discretization) when does a discrete Toeplitz determinant with Fisher- 
Hartwig singularities have the same asymptotics as a continuum one. Our 
result concerning this is that if there are N particles in the gas and M points 
in the discretization of the unit circle, and if N/M —)• 0 as —)• oo, then the 
discrete and continuum Toeplitz determinants with FH-singularities have 
the same asymptotics. Perhaps more interestingly, our analysis suggests 
that if N/M q > 0 as N ^ oo, the asymptotics of the discrete Toeplitz 
determinant will deviate from the continuum one - see the last section of 
this note for a discussion about this. Another interesting question (which 
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we don’t discuss any further) is how our results relate to discrete Riemann- 
Hilbert problems (see e.g. [2] for how these are used for studying asymptotics 
of discrete orthogonal polynomials). 

The outline of this note is the following: we start by describing our model 
and quickly review a construction for the relevant chaos measure. We then 
introduce our main result and the smaller ingredients this will consist of. 
Next we go on to review continuum and discrete Toeplitz determinants. We 
then prove our main tool for analyzing discrete Toeplitz determinants - a 
generalization of a result in [3], namely a description in terms a continuum 
Toeplitz determinant with the same symbol and a Fredholm determinant. 
Due to this representation, the relevant question then becomes analyzing the 
asymptotics of the Fredholm determinant which can be done by making use 
of the asymptotics of the (continuum) orthogonal polynomials with respect 
to a Fisher-Hartwig weight. The bulk of the note consists of straightforward, 
but tedious analysis of the asymptotics of the polynomials and the Fredholm 
determinant - this part relies heavily on results in [23, 24, 6]. Finally when 
we’ve proven the relevant asymptotics, we are able to prove our main re¬ 
sult (from the point of view of Gaussian multiplicative chaos it is Theorem 
9, but from the point of view of discrete Toeplitz determinants with FH- 
singularities it is Proposition 18). After the proofs, we discuss briefly the 
possibility of deviating from the continuum FH-asymptotics. 

Acknowledgements: The author wishes to thank Y.V. Fyodorov for 
discussions which have been extremely influential on the author’s view of 
log-gases and seminal for the work in this note. The author also wishes to 
thank A. Kupiainen and E. Saksman for discussions related to trace class 
operators, Hilbert-Schmidt operators, and Fredholm determinants. 

2. Model, main result, and outline of approach 

We will now review our model, introduce the relevant objects and notation 
to state our main result, and hnally discuss what are the relevant estimates 
we shall need. 

2.1. The model - a discrete log-gas and its characteristic polyno¬ 
mial. We begin by considering a discretization of the unit circle (M evenly 
spaced points on the unit circle) and then consider a probability measure 
on the A-fold product of this discrete unit circle with itself. The result¬ 
ing probability measure can be seen to be a discretization of the law of the 
eigenvalues of a Haar distributed N x N random unitary matrix. 

Definition 1. Let us denote by T, the unit eircle in the complex plane. Fix 
A'l G let • G —^ C, 

( 1 ) iMiz) = -I, 

and let 


(2) 


Vm = {z gC: '^m{z) = 0} C T. 
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We then fix some N G Z_|_ such that N < M and consider the following 
probability measure on 


(3) 


where 


Zn) 


1 

Zn,m 




(4) Zn^m = ^ \zi-Zj\^. 

We’ll denote or just K the expectation with respect to this measure 

(we’ll also drop the subscripts N,M from Fn,m when convenient). 

Remark 2. We’ll prove shortly that 


(5) 


Zm,m = N'.M^. 


Remark 3. The restriction that N < M is essential as Z^^m = 0 ifN > 
M (if N > M we can’t find N distinct points from Vm which contains 
M points). On the other hand, the other extreme would be when we let 
M ^ oo as N remains fixed. In this case we would simply get a log-gas 
on the unit circle and the law of the points {zi, zn) would be that of the 
eigenvalues of the Circular Unitary Ensemble (or Haar distributed random 
unitary matrices). 


The way we wish to describe the limiting properties of these probability 
measures as A^, M —oo is through the ’’characteristic polynomial” of the 
random points {zi,..., zn)- For simplicity, we focus on its absolute value (as 
the zeros give the locations of the points zt, this describes the global geom¬ 
etry just as well as the full characteristic polynomial). Also with the same 
effort, we can study (small enough) positive real powers of this quantity, 
so we add a further parameter to the quantity. Moreover, as we wish to 
prove convergence to a multiplicative chaos measure, we want to treat it as 
a measure so we make the following definitions: 

Definition 4. Let (3 > 0 and : C —)• [0, oo), 


N 

( 6 ) = Y[\w - Zjf. 

i=i 

Moreover, consider the following Radon measure on the unit circle: 


(7) 




dO 


We now describe what ^ will converge to as N, M 
way as well as the mode of convergence. 


oo in a suitable 
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2.2. Gaussian multiplicative chaos and convergence in distribution 
with respect to the weak topology of Radon measnres. The relevant 
limiting object is a Gaussian Multiplicative Chaos measure. We refer to [26] 
for a general review on the theory of such measures and to [4] for a concise 
proof for the existence (and uniqueness) of them. 

Definition 5. Let he i.i.d. standard complex Gaussians, L G Z_|_ 

and Xl : 27r) —)• M, 


( 8 ) XL{e)=-RX-Lzje»‘. 

For /? > 0, consider the Radon measures 


(9) 

„ e,hXL{e) ^0 

and 


(10) 

P^{d0) = lim p^^idO), 

L—)-oo 


where the limit is in the almost sure sense and with respect to the topology 
of weak convergence on the space of Radon measures on the unit circle, i.e. 
for every continuous / ; T —)• M 

r27r r27r 

(11) / f[F^)A{de) ^4“ / /(e*^)/(d0) 

Jo Jo 

almost surely. 

Remark 6. One can check that 

(12) E{XL{e)XL{e')) log - F^'\, 

so one can view Xl as a (constant multiple of a) regularization of the two- 
dimensional Gaussian Free Field restricted to the unit circle. If we write 

X for this field, then formally ^p{d6) = This expo¬ 

nential is purely formal as X is a generalized function and E(X(0)^) = oo. 
The proper definition of is precisely in terms of this limiting procedure. 
In [4, 19] it is proved that constructing pL^ through almost any other rea¬ 
sonable regularization of the field X yields a measure with the same law as 

Remark 7. A basic fact in the theory of Gaussian Multiplicative Ghaos 
is that indeed converges almost surely to a non-trivial random measure 
(which we call ) for 0 < /3 < 2. For /3 > 2, = 0. 

When proving convergence of the topology we’ll consider conver¬ 

gence in is still that of weak convergence of measures, but the convergence 
is no longer almost sure, but that of convergence in distribution. We refer 
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to section 4 of [21] for a general version of the following characterization of 
convergence in distribution to with respect to the weak topology of measures 
on the unit circle. 

Proposition 8. A sequence of random measures on the unit circle con¬ 
verges in distribution with respect to the topology of weak convergence of 
measures, to the measure v if and only if, as n —)■ oo, 

r27r r27r 

(13) / f{e^>n{de) A / fie^^Mdd), 

Jo Jo 

for every continuous /;¥—)• [0,oo). 

2.3. Main theorem. We are now in a position to state our main theorem. 
Theorem 9. Let N/M —)• 0 as —)• oo, and fJ G (0, y/2). Then as N ^ oo. 


(14) 4,M(rf^) ^ L^idO) 

in the sense of convergence in distribution with respect to the topology of 
weak convergence of measures. 

Remark 10. Note that we do not prove convergence for all values of 13. This 
is due to the fact that our approach relies crucially on calculating variances 
which are finite only for (3 G (0, \/2). It is an interesting open problem to 
extend such results to ft £ [\/2, 2) (or even further). 

Remark 11. ^s noted in Remark 3, one must have N < M to have a mean¬ 
ingful model, but this result still leaves open for example the cases where 
N/M —>■ g G (0,1] as N ^ oo. In fact, there are some hints in our approach 
that perhaps something else happens and the discreteness of the model (de¬ 
viation from the M —>• oo for fixed N - situation) becomes visible in the 
characteristic polynomial. For further discussion, see the last section of this 
note. 

2.4. Structure of the proof. We will now provide an outline for our proof 
of Theorem 9. Our argument is a rather standard probabilistic argument - 
we add a further level of approximation to by truncating the Fourier 

series of log ^(e*^), and we show that as N —)• oo, this converges in law to 
something that converges to as we remove the truncation. In addition to 
this, we prove that the variance of the error of this approximation tends to 
zero as we hrst let iV —)• oo and then remove the truncation. More precisely, 
let us make the following dehnitions: 

Definition 12. For L G Z_|_, let 


(15) 

and 


F\ 


'3 (pid\ _ 




(16) 
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The first (and easier) ingredient for the proof will be to show that ^ ^ 
converges to as —)• oo. More precisely 

Proposition 13. For any (3 > 0, L € Z_|_, and N, M such that M — N ^ oo 
as N ^ oo, 

(17) ^ 

as N ^ oo (with respect to the topology of weak convergence of measures on 
the unit circle). 

This will actually follow from convergence of certain linear statistics. 

Proposition 14. Let {Zj)^^ he i.i.d. standard complex Gaussians and for 
j e let 

N 

( 18 ) 

k=l 

Then for any fixed I G Z_|_, and any M,N such that M — N —)• oo as 
N —>■ oo, 

(19) {Zu Zi) A {Zi,V 2Z2, AlZi) 

as N ^ oo. 

Remark 15. Note that the conditions on (3, N, and M are mueh weaker 
here than in our main theorem. The faet that we don’t have any restric¬ 
tions on fd is due to the regularized field Xl being a non-singular object so 
exponentiation is trivial (the convergence of linear statistics is essentially 
a statement about convergence of the Fourier coefficients of log F^ m L 
those of the field {dX^). The faet that we have less restrietive eonditions on 
N and M is due to the faet that the charaeteristic polynomial seems to be 
more sensitive to the discreteness of the model than the linear statisties. 

The next ingredient for our proof will be estimating the variance of the 
approximation due to introducing L. 

Proposition 16. For any (d € (0, \/2), any continuous f :T ^ [0,oo), and 
N, M such that N/M 0 as N ^ oo, 


f(^"^)iA,M,LidO) - A,MidO))^ j =0- 

Expanding the square in the proposition above, we see that this boils 
down to precise, uniform estimates on (possibly mixed) moments of F^ ^ ^ 
and T)vm- relevant estimates are given in the following result. 

Proposition 17. Let (d G (0, \/2) and L G Z_|_. 

1) If M — N ^ oo as N ^ oo, 


C. WEBB 


(21) lim = 

N^oo ’ ’ 

uniformly in 6. 

2) If N/M —)■ 0 as N ^ oo, 


( 22 ) 


where G is the Barnes G function and o(l) is uniform in 6. 

3) If M — N ^ oo as N ^ oo, 

(23) E (<M,L(e*')<M,L(e*'')) = ^7V,l(0,0') X (1 + 0(1)), 

where o(l) is uniform in 9,9' and as N ^ oo, EMr(9,9') increases to 

4) If N/M 0 as N ^ oo, then 


^2 G (1 + I 


E{F^ 




id' 


N,M,L\ 


N,M 


1 . 

— ^ 4 2-^j=l j g 2 ^j=l 


■cos{j{e-9')) 


(24) 


X IV 4 


G{l + I 


G(l + /3) 


X (1 + o(l)), 


where o(l) is uniform in 9,9'. 

5) If N/M 0 as N ^ oo and there is a fixed 5 > 0 such that \9 — 9'\ > 5, 
then 


(25) 


G(l + /3) 


X (1 + 0(1)), 


where o(l) is uniform on \9 — 9'\ > <5. 

6) There is afixedtQ such that if N/M 0 as N ^ 00 , then for \9 — 9'\ < 
2to 


E [FN,Mie'^)FN,Mie'^')) = N^ 


2 G(1 + /3)2 

G(l + 2/3f 



(26) 


X e 


a2 2 sin 

-+log 


ie-e'\ 


\s-e>\ 


X (1 + 0(1)), 


where o(l) is uniform in \9 — 9'\ < 2io and a is a continuous function 
(depending only on ft - not on N, M, 9,9', or to ) on —with the following 
asymptotics: there exists a 5 > 0 such that 
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(27) = + 

as s —>• 0 along the negative imaginary axis and 


(28) a(s) = 0(|s|-') 
as s —>■ oo along the negative imaginary axis. 

The way these results are obtained is through a discrete Toeplitz deter¬ 
minant representation of moments of ^ and F^ We then show, 

mimicking an argument for non-singular Toeplitz determinants given in [3], 
that this discrete Toeplitz determinant can be represented in terms of a 
standard, or continuum, Toeplitz determinant with Fisher-Hartwig singu¬ 
larities as well as a certain Fredholm determinant. The asymptotics of the 
continuum Toeplitz determinants are well known (see [7, 8, 6] and references 
therein), so the problem then boils down to the analysis of the asymptotics 
of the Fredholm determinant and this is done using asymptotic proper¬ 
ties of orthogonal polynomials with respect to the Fisher-Hartwig weight 
on the unit circle (which have been worked out in or follow from results in 
[23, 24, 7, 8, 6]). One of the main results of this note is describing the regime 
for certain discrete Toeplitz determinants with FH-singularities where the 
classical Fisher-Hartwig formula still holds. We thus write the statement 
down explicitly (statements 2), 4), and 5) are essentially special cases of this 
result). 

Proposition 18. Let V be a Laurent polynomial of the form 

^ 1 . 

(29) V{z) = '^-{ajZ^ +ajZ 

3=0 

let be fixed distinet points on the unit cirele, and let jdj > 0 for 

j = 1,..., m. If we write for \z\ = 1, 


(30) f{z) = e^^^^\[\z-Wj\f^^ 

i=i 

and 


(31) 


r^_i(/) = det 





N-l 


p,q=0 


then as N ^ oo and M > N 
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T^_i(/) = e^^—+4 




(32) 


X n 

l<p<q<m 


Wp-Wq\ 2 


e 2 ^i=o( 2 ^j+ 2 %■ ) -f 

i=i 

2 

-) 

M J 


^ G 1 + ^ 

TT —^ 


1 + 0 


Remark 19. So in particular, if N/M —)■ 0 as —>■ oo, we recover the 
classical Fisher-Hartwig formula. If on the other hand N/M —^ G (0,1]; 'll 
is possible that this formula gets modified (for q small enough, our analysis 
indicates that it shouldn’t get modified too much - see our discussion at the 
end of the note). 


3. Toeplitz determinants 

We will now review some facts about Toeplitz determinants (discrete and 
continuum) as well as prove a connection between discrete and continuum 
Toeplitz determinants with Fisher-Hartwig singularities. 

3.1. Toeplitz determinants and orthogonal polynomials. We will now 
recall some basic facts about orthogonal polynomials on the unit circle and 
their relationship to Toeplitz determinants. 

Definition 20. Let / : T —>■ [0,oo) and let be the orthonormal 

polynomials with respect to the weight /(e*^)|^, i.e. (j)j{z) = XjZ^ +0{z^~^) 
with Xj 7^ 0 and 

r2ir _ jfl 

(33) <^i(e*')4(e*^)/(e*^)—= 

Remark 21. It is a well known fact that such orthonormal polynomials have 
a determinantal representation and from this determinantal representation, 
one can check (e.g. using the Heine-Szego identity) that for a non-negative 
weight which is not Lebesgue almost everywhere zero on the unit circle, the 
orthonormal polynomials exist and are unique. Moreover, one can check 
with this argument that for a real weight, Xj Is real. For more information 
see e.g. [7]. 

Let us now recall the definition of a continuum Toeplitz determinant. 
Definition 22. For a function / : T —)• C, we write 

(34) . 

\Jo 

Its connection to the orthonormal polynomials is given by the following 
result. 
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Proposition 23. Let / : T ^ [0, oo) be an function and let Xj be as in 
Definition 20. Then 

7V-1 

(35) TN-i{f)= n X-'. 

i=o 

We will also make use of the Christoffel-Darboux identity. 

Proposition 24. Let f be a weight on the unit circle such that the orthonor¬ 
mal polynomials exist and let denote the polynomial obtained by complex 
conjugating the coefficients of fij (i.e. 4>j{z) = fijifz)). Then for z,w G C, 
z,w 0, and z w 


(36) 


^ (t>j{w)fij{z ^^- 

i=o ^ 


Let us also recall the classical strong Szego theorem (see e.g. [31] for an 
extensive discussion related to the theorem). We consider it only in the 
generality we make use of. 


Theorem 25 (Szego). Let V be a Laurent polynomial of the form 


(37) Viz) = ^^{ajZ^ +ajZ ^). 

j=o 

Then 

(38) + i + „(!)) 

as N ^ OO. 


Finally let us recall the asymptotics of Toeplitz determinants with Fisher- 
Hartwig singularities. The results in the form we state them can be found in 
[33, 8, 6], though for a more comprehensive picture of the problem, see also 
e.g. [7] and references therein. For the uniformity in the first statement, see 
e.g. [8] where it is mentioned (it is a consequence of the uniform estimates of 
the jump matrices in the ”hnal transformation” and the uniform expansions 
of the relevant local parametrices). 

Theorem 26 (Widom; Deift, Its, and Krasovsky; Claeys and Krasovsky,...). 
1) Let V be a Laurent polynomial of the form 

^ 1 . 

(39) Viz) ='^-{ajZ^ +ajZ ^). 

Moreover, let iwfiJLi be distinct points on the unit circle such that for 
some fixed d > 0, \wi — Wj\ >6 for i j, and > 0 for j = 1,..., m. Then 
for 
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(40) 


j=i 


Tn-iU) = H ^(«*-‘+^%-‘)ArE: 

j = l 


j' = l 4 


(41) 


- n 

l<p<q<m 


_ 

Wp - Wg\ 2 


^ g(i + 4 


n 


-(i+°(i)), 


^(1 + /^i) 

where o(l) is uniform on \wi — Wj\ > h and G is the Barnes G-function. 
2) For t > 0, \z\ = 1, and /3 > 0, let 


(42) ft{z) = \z- FY\z - e-^Y- 

There is a fixed to such that for t < to 


(43) 


TN-iUt) = 


G{1 + fif i 
G{l + 2fif 



X e 2 t X (1 


0 ( 1 )), 



ds 


where o(l) is uniform in t < to and a is a continuous function (depending 
only on fi - not on N,t, or to) on —iK_(_ with the following asymptotics: 
there exists a 6 > 0 such that 


(44) ^^s) = ^fi^ + 0{\sf) 

as s —)• 0 along the negative imaginary axis and 


(45) <j(») = 0(|8|-‘) 

as s —>■ oo along the negative imaginary axis. 

3.2. Discrete Toeplitz determinants. A discrete Toeplitz determinant 
is one where the integrals appearing in the entries of TN-iif) are replaced 
by an approximating Riemann-sum. We now argue how these arise from 
expectations with respect to the measure 

Proposition 27. Let / : T —)• C he an arbitrary function. Then 


N 


(46) n = TN-i{f) ■■= det I - 


i=i 


zGT>m 


N-l 


j,k=0 


Proof. The proof is essentially the same as in the continuum case, but we 
write it down for the convenience of the reader. By definition 
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(47) 


N N 

j=l ’ ^ zi,..,2jveX>M ^<i<j<N k=l 


We then make use of the fact that the hrst product is simply the square 
of the absolute value of the Vandermonde determinant: multiplying out the 
two determinants and using multilinearity in the first row we hnd 


n 1 ^* “ 

1<2<J'<W 


(48) 


1 

1 

1 


1 



Zi 

^2 

ZN 


1 

^ . 



JV-1 . 
^2 

Zn 


1 

w ■ 



1 ^2 * * * '^N ^ 

N ••• 

Ef=i^. N ... E,=iE^-^ 


Af-l N-1 v^A-1 N-2 


-yrN-l 
2^j=l ^3 

—N-2 

2^j=i 


N 


N 


E 

1=1 


^3=1 V 
1 
V 


^3=1 V 

Ej=i 

N 


Z; 


N-l 

I 


s;^N—l N—2 
^i=i ^3 


Now in the sum in the second row, the Itli term in the sum is proportional 
to the hrst column (the factor of proportionality being ^) thus - making use 
of multilinearity - this term in the sum yields a determinant equal to zero 
and we see that 


(49) 




1 

ZI 2 

N^N — A-1 

l^i=i V 

n \zi-zj\‘^= 

N 

E E 

Vi 

1 

y^A —N-2 
■ ^j=l V 

l<i<_7<A^ 

/l=l l<l2^N,l2^ll 


N-2 

V 2 

N 

Continuing in this 

manner, we hnd 





(50) 


n 1 ^* “ 


E 

crSSjv 


1 Zai2) 

^cr(l) f 


N-l N-2 
^a(l) ^a{2) 


Za{N) 

Z(t[N) 


■N-l 

■N-2 


1 


Now as n^i fi^j) is invariant under permutations of the indexes j, we 
see by multilinearity that 
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N 


En/( 


Zi = 


J=1 


(51) 


A^! 

Zn,m 


E 




1 ^^2 

Zl 1 

N-1 N-2 


Nl 

Zn,m 


det 


iZG'DM 


N-1 


j,k=0 


Then setting f{z) = 1 for all z G C we see that 


ZN 

W 


■N-1 

■N-2 


N 


n 

j=i 


(52) 


1 = 


A^! 

Zn,m 


^M-l 


n-1 


,-{3-k)2-Ki- 




\ l—O / j,fc=0 

For j ^ k, the sum in the above determinant is zero and for j = A: it is 
Ad so we find 


(53) 

which implies that 


Zmm = N\M^, 


N 

(54) lEiV,Mn/(^i) = ^^-i(/)- 


Remark 28. Note that as we are interested in moments of ^ ^ and 

^N M’ symbols f that we are interested in are special cases of the Fisher- 
Id artwig symbol: / : T —)• C, 


k 

(55) f{z) = e^^^'>Y{\z-Wj\^i, 

j=i 

where V is a smooth enough function, are fixed distinct points on 

the unit circle, and fij >0. In our case, V is a Laurent polynomial, either 
k = 1 or k = 2, and either (3j = 0 or fij = /3. ITe will also need the 
situation where /3j = 0 for all j. When it is not restrictive, we will consider 
the general symbol (55). More generally, one could also consider the phase 
of the characteristic polynomial (causing the symbol to have a jump), one 
could consider more general potentials V, or complex /3j. 

3.3. Connection between discrete and continuum Toeplitz deter¬ 
minants with Fisher-Hartwig symbols. We now present our main tool 
for analyzing discrete Toeplitz determinants with Fisher-Hartwig singular¬ 
ities - namely a generalization of an argument of [3] where it was noted 
that the discrete Toeplitz determinants with no singularities can be writ¬ 
ten as a product of a continuum one with essentially the same symbol, and 
a Fredholm determinant related to the continuum orthogonal polynomials 
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and a function related to the discretization. Compared to their approach, 
our symbol does not have an analytic continuation into a neighborhood of 
the unit circle, but it does have a continuation with branch cuts along the 
rays going through the singularities. Our approach is to modify the argu¬ 
ment of [3] to this case. We begin by recalling this continuation (for more 
information, see [24, 7]). 

Lemma 29. Let V be a Laurent polynomial of the form 

p 

(56) V{z) = 

j=-p 

let be distinct points on the unit circle, and for \z\ = 1, let 

k 

(57) f{z) = e^^^^ll\z-w,fL 

i=i 

with Pj > 0. There is an analytic continuation of f into C\{Lij^iWj x (0, oo)) 
and it can be written as 


(58) 


/(^) 


A(/, z) 
De{f,zy 


where Di{f,z) is the analytic continuation of 


(59) z i-A JJ(1 - zw^y^ 

1=1 

from {\z\ < 1} into the domain C\ x [l,oo)) (for \z\ < 1, the root 

is according to the principal branch). 

Similarly, De{f,z) is the analytic continuation of 


(60) Z 

1=1 

from {jzj > 1} into the domain C\ x [0,1]) (again in \z\ > 1 one 

takes the principal branch). 

Proof. For \z\ y 1, consider the function 


(61) 


D{f,z) = exp 


r-27r 


logfie 


— z 


dO 


For \z\ < 1 we have 


1 

47r 
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log D{f,z) = 


1 /*27r / P ^ 

■^0 \j=-p 1=1 


e*® - wi\ 


(62) X (1 + e-^^z) e-^‘i^z'ide. 

q=0 

Noting that for j S Z \ {0} 


(63) 


2vr Jq 


log |e*® - w\e^^Ue = - —w^ 


m 


and log |e*® — i(;|d0 = 0 so that we get 


(64) 


or 


iogo(/..) = |+E%-'-ETEl 


k ^ oo ^ 

Pi i 


i=i 


2 ^ 


i=l j=l 


(65) 

= e^+E^.i a^P+Eti ^ iog(i-^«<r') 


n-z^ 

= P 2 ^ ^j=l 


n(i 


— zw 


-1 


i=l 


2 


where the branch of log(l — zw^^) is chosen by the series expansion. This 
function is analytic in \z\ < 1 and extends analytically into C \ x 

[l,oo)). We will denote this continuation by Di(f,z). 


Similarly for \z\ > 1 we have 


1 p 27 V I P ^ 

log D{f, z) = -— ^ + X] 

^ -'0 \j=-p 1=1 

OO 

(66) X + 

q=0 

p k a OO 

= -|i-E«-w' + E§Er-'P 

3=1 1=1 3=1 

or 

D{f,z) - 

(67) = ]^(i _ 

1 = 1 

where the branch is given by the series expansion. This is again analytic 
in \z\ > 1 and continues analytically into C \ {Uj^iWj x [0,1]). We call the 
analytic continuation De{f,z). Then 
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(68) 


f{z) 


Deif,z) 


is analytic in C \ x [0, oo)) and agrees with the original symbol on 

the unit circle. 

□ 


We are now ready to adapt the proof of [3] to prove the connection between 
the discrete and continuum Toeplitz determinants with Fisher-Hartwig sin¬ 
gularities. We restrict to the case where V is real on the unit circle simply 
for the reason that this is the case we need. 

Proposition 30. Let 1/ : C \ {0} C be a Laurent polynomial of the form 


p 

(69) V{z) ='^^{ajZ^ + ajZ~^), 

j=0 

{wi, ...,Wk) be distinct points on the unit circle such that 0 < arg(t(;i) < 
arg(u; 2 ) < ... < a,\:g{wk) < ^tt, and for \z\ = 1, 


k 

(70) f{z) = e^^^^\[\z-w,\f^L 

j=i 

where (3j > 0. We also write f{z) for the continuation of this function 
described in Lemma 29. 

Let e G (0,1), and let us write Fj for the contour consisting of the segments 
{\ — e,l + e)xwj, (1 — e, l-|-e) X rcj+i and the (shorter) circular arcs of radius 
1 ± e joining (1 ± e)wj and (1 ± e)wj+i. For j = k, we understand Wk+i as 
wi- We also understand the contour to be oriented in the counter-clockwise 
direction. 

Moreover, let 


(71) 


v(z) 


l_zM ) 
Z-M 

l-Z-M ’ 


z\ <1 
z| > 1 ■ 


Finally we denote by a the measure on UjFj defined by 


(72) 


a{dz) 



z\ = l±e 
z\ G (1 — e, 1 -|- e) 


and by (Bj^iL‘^{Tj , a) the space of functions h = {hi ,..., hk), where hj : Fj —>• 
C are square integrable with respect to a. 


Then we have 


(73) 


TN-i{f) = TN-i{f) det(l 
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where K is a trace class integral operator on (BjL'^{Tj,a) with kernel (acting 
on each component of the vector) 


(74) K{z,w) = I{z)y/v{z)f{z)y/v{w)f{w)KcD{z,w), 

where the roots are according to the principal branch, when acting on the 
jth component the values of f{z) and f{w) are calculated as limits from the 
inside ofTj, 


(75) 

and 


I{z) 


1, |z| = 1 zb e 

-i, \z\ G (1 - e, 1 + e) ’ 


Kcd{z,w) = -—-. 

Z 

Remark 31. Note that we can think of an element h G as 

defining a two-valued function on - or more precisely, it’s single valued 
on (1 zb e)T and two-valued on (1 — e, 1 + €)wj. 

Remark 32. Keeping in mind the previous remark, an important difference 
between our case and the non-singular case considered in [3] is the different 
form of the integration contour. Instead of just circles of radius 1 zb e, we 
also have the segments joining these and as there is a branch cut along 
each such segment, the oppositely oriented integrals don’t cancel and we will 
need slightly more eomplicated asymptotics of the orthogonal polynomials to 
control our Fredholm determinant. If = 0 for all j, there are no branch 
cuts, the oppositely oriented integrals along the rays through Wj cancel, and 
we recover the result of [3]. 

Proof. The proof is essentially a simple modification of that in [3] to take 
into account the branch cuts of /. To do this, we being by noting that as 
f{wj) = 0, we have 



(77) 


TN-i{f) = det 



z£VM\{wi,...,Wk} 


-{l-m) 


fiz) 


We then note that by Cauchy’s integral formula, 


N-l 


l,m=0 


(78) 


r^_i(/) = det 



. 7 m (^) dz 


N-l 


/,rn=0 


where on Tj, / is defined by a limit from the inside of the contour. 


One might need to take some care here as the contour Tj goes through 
the points wj, wj+i so if wj G Vm, the integrand might have a singularity as 
lM{wj) = 0 in this case (so at most, a simple pole). But as also f{wj) = 0 
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(or more precisely, f{z) = 0{\z — Wj\^^) as z Wj) we see that this is an 
integrable singularity so one can indeed apply Cauchy’s integral formula. 

As determinants are unchanged by elementary row and column opera¬ 
tions, this implies that 


1 


Tn-iU) = ^N-l 2 
llj=0 X.j 


(79) 




, JT, 


A 7m(^) 

\ N-1 
dz \ 

M'Jm{z 

) 27ri 1 

/ /,m=0 
\ N-1 


dz 1 

M-fMiz) 

to 1 

. 1 

II 

O 


where we made use of Proposition 23. Now using the orthonormality con¬ 
ditions, we write 


(80) 


J 

= ^l,m - f Mz~^)(j)miz)f{z) 

J\z\=l 


dz 

27riz 


+j:/ 


As 

(81) 


Z^ct)i{z ^)(j)m{z)f{z) 


1 


is analytic inside of Tj , we can deform the part of T that is inside of Pj into 
the curve {z G Pj : \z\ > 1} with opposite orientation (by Cauchy’s integral 
theorem). Thus we find 




(82) 


dl,m + 


r,n{|2|<i} 


4>i{z )(pmiz)fiz) 


^7m(^) dz 
M'^m{z) 2'Kiz 


J Jrjn{\z\>i} \M'^m{z) ) 27 rz 

= + M^~^)^rniz)f{z)v{z)- 


2'niz 

3 - 

Now we need to write the determinant as a Fredholm determinant. Let 
us introduce the following notation: let 


(83) 
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(84) {A{hi,...,hk))i = '^J^ (j)iiz ^)^/v{z)f{z)hjiz)^— 
and 

(85) B : f ({0,iV - 1}) ^ {Tj, a ), 

N-l 

(86) S(ao, ...,aAr_i) = ^ 

k=0 

where both of the roots are according to the principal branch, and the jth 
component of this function is obtained by restricting to Fj and defining / 
as a limit from the inside of Fj. 

Note that by Cauchy-Schwarz, since Pj > 0 (or more precisely since 
/(2) = O(|2-t0j|ft)), (A% is indeed finite for any h £ jic) even 

if V has a simple pole at some wj (which occurs when wj £ T)m)- For the 
same reason, B really maps into ©jL^(Fj,u). 

Thus we have 

( 87 ) {AB)i^k = '^ [ '4>i{z~^)(l^k{z)v{z)f{z)^^ 

and 


(88) T^-i(/) = Tn-iU) det(/ + AB). 

We note that as K = BA is an operator on and 


r - . _ dz 

(89) Kh = Y.Y. ^)^/v{z)f{z)hj{z)^^f<j)k, 

k=0 j 

(where the jih. component is obtained by a limit from the inside of Fj) so 
we see that K a finite rank operator, and in particular, trace class. 

By Sylvester’s determinant identity, det(/ + AB) = det(I + BA) (/ de¬ 
notes the identity in the relevant space and the determinant is a Fredholm 
determinant). 

Making use of the Christoffel-Darboux identity (Proposition 24), we find 
that the kernel of K {{Kh){w) = ^ . L, K{z,w)hj{z)a{dz)) is 

'' 3 


(90) 

K{z, w) = I{z)^/v{z)f{z)^/v{w)f{w) 


(pN{z)4>Niw ) - 4>Niz A4 >n{w) 


1 - 


□ 
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Remark 33. A similar result and argument also holds for a log-gas on a 
discrete subset of the real line, where one ends up with a discrete Hankel 
determinant, which can be related in a similar way then to a continuum 
Hankel determinant. 

4. Asymptotics of orthogonal polynomials 

We will now review the basic results of the asymptotics of the orthogonal 
polynomials with respect to the Fisher-Hartwig weight on the unit circle 
(as well as the non-singular case) from [23, 24, 6]. We will consider the 
asymptotics on the contour relevant to us. We will start with the non¬ 
singular case studied in [23]. Next we will consider the singular case with 
the restriction that the distance between singularities is bounded away from 
zero (studied e.g. in [24]). Finally we will consider the situation where the 
singularities can merge, i.e. their distance can go to zero as N ^ oo. This 
case has been studied in [6]. As in some cases the exact estimates are not 
written down explicitly in these articles, we will review a small part of the 
relevant proof, to argue why the estimates we require hold. 

4.1. The non-singular case. We will hrst review results from [23] which 
we will use to analyze the asymptotics of the Fredholm determinant in the 
non-singular case. This has already been essentially done in [3], but as in 
the proof of Proposition 16, we wish to integrate the resulting (discrete) 
Toeplitz determinant, we’ll need uniform asymptotics, so we will take some 
care. 

Theorem 34 (Martinez-Finkelshtein, McLaughlin, and Saff). Consider a 
Laurent polynomial V of the form 

(91) V{z) = +ajz-^), 

1=0 

let {(j)j)'^Q, with (j)j{z) = 4- 0{z^~^), be the orthonormal polynomials 

with respect to the weight e^^^^ on the unit circle (as > 0 on T, the 

polynomials exist). Also fix any R > 0 and e G (0,1). Then there exists a 
c > 0 such that for \z\ = 1 — e, 

(92) (fwiz), = 0(1)6"''^ 

as N ^ oo. Moreover, the 0(1) term is uniform on (1 — e)T and on the set 


(93) Ar = < max | a,- [ < R> . 

[i<l<p J 

For \z\ = 1 4- e, one has 

(94) <PNiz)=0{l)il + e)^ and </.)v(^) = 0(l)iV(l + e)^, 
where again the 0(1) terms are uniform on (1 -|- e)T and Ar. 
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Remark 35. The same results hold for the conjugate polynomials 4>j{z) = 

Proof. This is essentially an application of Theorem 1 and Corollary 2 in 
[23]. We will not go into great detail here, but try to point out where in 
[23] the reader can convince themselves about the validity of the required 
estimates. The only things that aren’t immediately clear from the results in 
[23] are the asymptotics of the derivatives and the uniformity on Aji. 

For the derivatives, we note that (in the notation of [23]) the derivative of 
has similar bounds as those in equation (24) of [23]. Then by uniform 
convergence, one can differentiate £n and term-wise, which gives the 
desired estimate for the derivatives. 

To have uniformity on A^, the bounds (24) in [23] (and corresponding 
ones for the derivatives) imply that it is enough to have uniform bounds on 
the term A defined in (20) of [23]. This in turn can be checked from the 
relevant definitions. Checking that the convergence of xn is uniform on Ajv 
is similar. 

□ 

4.2. The singular case, when the distance between singularities 
is bounded away from zero. Again the relevant asymptotics have been 
studied in [24], but the uniformity in V or the estimates for the derivatives 
are not stressed. Due to this, we’ll again sketch some of the arguments for 
checking these facts. 

Theorem 36 (Martfnez-Finkelshtein, McLaughlin, and Saff). Let F : C\ 
{0} —)• C 6e a Laurent polynomial of the form 

(95) F(z) = ^i(ajF To^z"^), 

3=0 

{wi, ...,Wk) be distinct points on the unit circle such that 0 < arg(t(;i) < 
arg(u; 2 ) < ... < aig{wk) < 27r, and assume that there exists some (smalle 
enough) fixed <5 > 0 such that \wi — Wj\ >6 for i j. For \z\ = 1, let 

k 

(96) f^,) = enAYl\z-w,fL 

i=i 

where fij > 0. We also write f{z) for the continuation of this function 
described in Lemma 29. Finally let {4>j)j be the orthonormal polynomials 
with respect to the weight f on the unit circle (as f is non-negative and 
almost everywhere non-zero these exist) and write (fjiz) = 0{z^~^). 

We then have 

1) For any fixed compact subset of the unit disk, 

()n{z) = 0{N-^) and cj'^iz) = 0{N-^) 
uniformly on the given compact set, as N ^ oo. 


(97) 
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2) For any fixed compact subset o/{|z| > 1}, 


(98) finiz) = 0{1) and fi'^iz) = Nz^O{l) 

uniformly on the given compact set, as N —>■ oo. 

3) On the interval [1 — 5/6,1 + 5/6] x Wj, 


(99) (t>N{z) 


and 


0 ( 1 ) 

0 ( 1 )- 


y/ fO) 


m-\zm 


VfO) ’ 

fo(l)iV^, 

\om-\z\ 


|1 1^11 < N 

|1 “ 1^11 > jf 


|1 - Ull < 


N 


\o\ > N 


( 100 ) = 


0(l)iV^(iV|l-|z||)- 

fo(l)iV^+\ 

\o(i)|i-N|-^- 


|i kll < N 

|1 “ 1^11 > F 


|i - Ull < 


N 


|1 kll > N 


where the 0(1) terms are uniform on the interval and Xj\Jf{z) is 
understood to be the limit of 1/ y/f(ze^ as e —>■ 0 (and the root is 
according to the principal branch). 

Moreover, if Ar is defined as in (93), the above estimates are uniform on 
An and they are uniform on Bs = {\wi — Wj\ >5 for i / j}. 


Remark 37. Again corresponding results hold also for the conjugate poly¬ 
nomials 4>j. 

Remark 38. We have two different representation for the asymptotics of 
the polynomials on the rays as in some cases we ’ll make use of the fact that 
there is the factor fff{z) appearing in the kernel of K so this will cancel 
with the factor o/ l/i/ f{z) when analyzing the asymptotics of the kernel. 


Proof. While the asymptotics of the polynomials follow from Theorem 1, 
Theorem 2, and Theorem 3 in [24], the asymptotics of the derivatives are 
not written out explicitly in [24] and the uniformity in V is not stressed. So 
let us review parts of their argument. 


Using (by now standard) Riemann-Hilbert arguments, it is proven in [24] 
that in a compact subset of {|z| < 1}, if we write fijsf = xn^n, then 
4>7 v(^;) = Yii{z) and Xw-i = —^ 21 ( 0 ) (our normalization of the weight 
differs from theirs by a factor of 27r), where T is a 2 x 2 matrix valued 
function which is analytic in C \ T and in the compact subset it can be 
written as 


(101) 


Y{z) = S{z)N{z) 
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where S{z) is a 2 x 2 matrix valued analytic function in the given compact 
set, and satisfies 

(102) S{z) = I + 0{N-^) and S'{z) = OiN-^) 

uniformly on that compact set. Here I is the 2x2 identity matrix and 
the estimates are entry-wise (a convention we will follow through this note 
unless otherwise stated). Moreover, from the Neumann-series representation 
of S (see Proposition 4 in [24]) it follows that these terms are uniform on 
An and Bs- More precisely, this involves obtaining uniform estimates for 
the jump matrix J 5 which in turn boils down to uniform estimates on the 
scattering function S. These can be checked from its definition. 

Inside the unit disk one has 

(103) Niz) = {D,{f,z)D,if,0)r , 

where 


(104) 

From this, we find 



(105) ^Niz) = z)Di{f, 0)r^Suiz) = 0{N-^) 

and 


(106) = 0{N-^) 

uniformly in the compact set, uniformly on An, and uniformly on Bs- 
Moreover, 

(107) x%-i = 522(0)A(/,0)-2 ^ A(/,0)-2 

uniformly in An and Bs ss N ^ 00 (so in particular, it is uniformly 
bounded). 

In a compact subset of {\z\ > 1}, one has the representation 

(108) Y{z) = S{z)N{z)z^^^, 
where again 

(109) S{z) = I + 0{N-^) and S'{z)=0{N-^) 
uniformly in the compact set and uniformly on An and Bs- For \z\ > 1, 

( 110 ) N{z) = {D,if,0)Deif,z)r- 

Thus 
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(111) = 5ii(z) A(/, 0)Deif, z)z^ = 0{l)z^ 

and 

(112) <l>'^{z) = 0{l)Nz^ 
uniformly on the compact set and uniformly on An and Bs- 

In neighborhoods of the points wj, the expression for Y(z) is slightly more 
complicated. In particular, the functions in terms of which Y is written 
have branch cuts (some of which cancel when taking the product). As 
‘hw nevertheless is analytic, we can calculate ^^(- 2 ) on the interval 2 : € 
[1 —(5/6, l + (5/6] X Wj as a limit of <l>jv(^0 as 2 ' ^ 2 ; in any way we wish. The 
way we shall calculate it is by writing z' = \z\wje^'' and then let e —?> 0^. For 
simplicity, we’ll focus on the | 2 ;| < 1 case. The |2;1 > 1 case can be analyzed 
in a similar manner. 

One still has ^Niz) = 111 ( 2 ;), but now, for z' = \z\wje'''^ for e > 0 small 
enough, T is written as 


(113) Y{z') = S{z')E{z')^ ( -i— log ^ 






0-3 


VK- 


-izr 


For such z', T is of the form 



Here is a modified Bessel function of the first kind and we denote by * 
the second column as it will be unimportant to us. 

Moreover E is analytic in the given neighborhood of Wj and 

(115) £( 2 ) = (Q).*))"’T (_\ i) 

where in the given neighborhood of wj, Q is a non-vanishing analytic func¬ 
tion which is independent of N and can be bounded uniformly on An and 
Bs (as can its derivatives). 

Finally S is as before, i.e. it is analytic in this neighborhood of wj, and 
it satisfies 


(116) S{z) = I + 0{N-^) and S'{z) = 0{N-^) 

uniformly in the neighborhood and uniformly on An and Bs- 
Writing out the definitions, we find that for | 2 | < 1 
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( 117 ) 


^n{z') = iQ{z')Sn{z') - QizT^S^2{z')) 


N . 


(118) xJ-i—log — 1/3^. , log— +7/3, 1 


Wi \ -t+2 


N . 




^3 _ 

2 2 


To simplify notation a bit, we note that 




(119) ^ \/f {Qi^')Sii{z') - Q{z')-^Si2{z')) e'T- 

is 0(1) (uniformly on Ar and Bs) as is its derivative. We thus see that if 
we write 

N z' 

(120) c. = -i7log—, 

2 Wj 

we have 


(121) ci,^(/) = 0(l)^^e-^^v^ (Ce) + 7/3, Aq) , 

Af{z') V ^+2 ^-2 / 

where the 0(1) term and its derivative are uniform on and We then 
take the e —>■ 0 limit of this and find (for C, = — A^(log \z\)/2 > 0) 

(122) ^ n { z ) = 0(l)-^e-«yC (^^+1 (0 + (O) , 

where the expression l/\/ f{z) is dehned precisely as this limit. 

We now make use of the asymptotics of 7,,: for 0 < x < 1, 


(123) 

and for X > 1, 

(124) y/xe~^ I^{x) = 1 + 0(x“^). 

Thus if — A^log \z\ < 2, (i.e. C S (0,1)) 

(125) <h^(z) = 0(l)-^(-iVlog|2|)^ = 0(i)^^(Ar|i _ |^||)^. 

VfA) VfA) 

Whereas for — A^log \z\ > 2, we have 


^^{z)=0{l)-^ 

v/(+ 


( 126 ) 
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As both of the expressions agree (up to a uniform 0(1) factor) when 
|1 — \z\\ = the difference between the conditions —A^logj^l < 2 

and say A^|l — |z|| < 1 is immaterial so we have proven the claim about the 
asymptotics of the polynomials. Let us now consider the derivatives. 

For this, we use the fact that the 0(1) terms had uniformly 0(1) derivatives 
and that Il(x) = ^(Ij^+i(x) + /y_i(x)). We also note that for ( = log 


(127) 


d_ _ ^ld_ 

dz 2 z dC 


If we write 1(0 = IPj j (C) + IPj i (Oi we conclude that (with a similar 
~ 2~^2 ~2 2 

limiting argument) for \z\ < 1 on our interval 


<l>0(z) = O(l)-^^e-‘^ 

(128) X (o(i)i(C) - l^KC) - + Li(c) - "i'(c)) . 

We then note that (uniformly) 

(129) Q^ = 0(l) + ^^. 

f(z) z - Wj 

Putting things together (and using the small x asymptotics of ly) we see 
that for Q G (0,1), (i.e. roughly A^|l — |z|| < 1) 


(130) 


^'n0) = O(1) 




--VCe ^ 


ok 


N - 


A- 


z-WjY(^ + l) 



1 

2 


3 ^3 3 

8^ r(§ + \) Uy 4^ r(§ -1) Uy 

where all of the estimates are uniform. We then note that 


m + \) 


ii_3 
2 2 


Pj C 

4 z — w-i 


= N- 




(131) 

8r(f+ 4) 

and that for x G (1 — (5/6,1), 


8z 4z V 2 2j J 



log X 1 
X — 1 X 


(132) 


0(x- 1). 
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Thus 


(133) 


r(| + i) 


--/ /tj c 

N 

N 

1 

V 4 Z-Wj 

8z 

4z 

\2 2 

we conclude that for 

Ar|l - 

- |2:|| 

< 1, 


= o c 


(134) 4>(v(^) = 0(l)^^iVC^ = 0(l)^^iV(iV|l - \z\\fi- 

V/(^) vfi^) 

For A^|l — l^ll > 1, we have by (128) and (124) 

( 135 ) 

l-'ivW =0(l)^i^ (o^T + 

\/fiz)\\\z\-i-\ ) 'J f\Z) 

uniformly on and B^. 

□ 

4.3. Merging singularities. The analysis in the case of the merging singu¬ 
larities becomes very involved quite rapidly, so we will focus on the simplest 
case which is sufficient for us - namely we consider F = 0, A: = 2, and 
/3i = /32 = /3 > 0. In [6], a more general case (still with two singularities) is 
considered, but we restrict to this case. Let us now review the implications 
of the analysis in [ 6 ] for the asymptotics of the orthonormal polynomials 
with respect to such a weight. 

Theorem 39 (Claeys and Krasovsky). Fix a small to > 0. Then for t G 
(0, to) and (3 > 0, let ft : T ^ M, 


(136) = 

Let {4’j)^0’ 4>j{z) = XjZ^ + 0{z^~^) be the orthonormal polynomials with 
respect to the weight ft on the unit circle. We then have the following as¬ 
ymptotic behavior. 

For any fixed compact subset of the unit disk, as N ^ oo 


(137) cj)Miz),(^%{z) = 0{l)N-^ 

uniformly in z and 0 < t < to. 

For any fixed compact subset o/{| 2 ;| > 1}, 


(138) 4>n{z) = z^ 0{1), and 4>'^{z) = Nz^O{l) 

uniformly in z and 0 < t < to- 

For any fixed 7 > 0 and z G [1 — 7 ,1 + 7 ] x the asymptotics depend 
on the behavior oft. For certain fixed c > 0 small enough and C > 0 large 
enough, we have the following behavior: 
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For Nt < C 


(139) 

4>Niz) 


and 



1 

Vfdz) 

1 


iV^|l 


1^11 2 max(|l 


1 ii 2 . 


0{l)Nf^, A^|l-|z||<l 

0{l)\l-\z\\-P, iV|l-|z||>l 


(140) 

( 0(1) } — |z|| 2 max(|l — jzll 2 ,^ 2 ), 

</,(y(z) = J V/‘(-) 

_f0(l)A^^+\ A^|l-|z||<l 

“ |0(l)iV|l - |2||-^, iV|l-|z||>l 

uniformly in z and 0 < t < cN~^. 


For cN-^ <t<CN 


(141) 


4>n{z) = 


and 


^-1 



-|z||)f, iV|l- 


A^|l ■ 

0(l)iV^, 

iV|l - |z|| < 1 

0(l)|l-|z||-^ 

A^|l - |z|| > 1 


r0(l)iAr(iV|l_ 1^11)2, N\l 
(142) (p'j^{z)={ 

_f0(l)iV^+\ iV|l-|2||<l 

“ |0(l)iV|l - iV|l-|2||>l 

uniformly in z and cN~^ < t < CN~^. 

t > CN-^ 


(143) (j)N{z) 


0(l)^(iV|l-|^||)§, iV|l-|zl| 
0(1)^, N\l-\z\\ 

0(l)(iVt-i)i 

0(1)|1 — \z\\~i min(|l — \z\\~i 


A^|l - |z|| < 1 
A^|l - |z|| > 1 


iV|l - |z|| < 1 

A^|l - |z|l > 1 


< 1 
> 1 


z|| < 1 

z|| > 1 


1 

1 

A^|l - lz|| < 1 
A^|l - |zl| > 1 
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and 


(144) 


(pNiz) 



iV|l- I 2 II < 1 


|1 - |z|| G {N- 


|1 — jzll > t 


t) 


< 0 ( 1)|1 — \z\\~2~^t~2^ 


iV|l - |z|| < 1 
|i - |z|l e 
|i — l^ll > t 


uniformly in z and CN ^ < t < to. 


Remark 40. The orthogonal polynomials with respect to the weight \z — 
wp\z — W 2 \^ with largwi — aigW 2 \ = 2t are obtained from the ones with 
weight ft by simply rotating the argument - in partieular, the same asymp¬ 
totics hold near Wi. 


Proof. In [ 6 ], it is not the asymptotics of the orthogonal polynomials but 
the Toeplitz determinant that was the main focus, so to check the required 
estimates requires a fair amount of work from us. 


The proof for the compact subset of the unit disk is roughly the same 
as in the case where the distance of the singularities is bounded away from 
zero. In [ 6 ], it is proven that there exists a function T (corresponding to S in 
the previous case) which is analytic in C apart from a jump contour which 
consists of the boundary of a fixed disk U containing the points for 
t < to (but not intersecting the given compact subset) and part of a lens on 
the unit circle going from e** to in the counter-clockwise direction - also 
not intersecting the compact subset (similar to the case when the distance of 
the singularities is bounded, except here we only have a single lens instead 
of one going also from to e** in the counter-clockwise direction. This 
function satishes T{z) = I + 0{N~^) and T'{z) = 0{N~^) uniformly off of 
the jump contour. 

It follows from the analysis in [ 6 ] that one then has for {|z| < 1 : 2 
(here we write r = Di{f,0)~^) 


(145) 


= Mf) = ( T(.) 


XN 


Diifuz] 


ca 


0 

-1 


11 


and xn-1 is given by in terms of the (2,1) entry evaluated at zero. From 
this one hnds the claim for the compact subset of the unit disk. 


In the case of the compact subset of {|z| > 1}, the argument is again 
analogous to the situation where the distance between the singularities is 
bounded. 


The part for z £ [1 — 7 ,1 4 - 7 ] x is more complicated. We now choose 
U so that it contains these segments. For simplicity, let us consider the 
situation where z = (1 — r) x e®* for r > 0. The other cases can be treated 
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in a similar manner. Moreover, we calculate things for such a z by taking 
the e —7- 0^ limit of (1 — r) x In this region, we can write 


(146) c^^(z) = (T(^)P(z))i,, 

where T again satisfies the same conditions and 


P{z) 


^ Di{ft,z)De{ft,z) 


(147) X ^ 




4- 


log z, —2iNt 


and 41 is defined by the following Riemann-Hilbert problem (in [6], it is 
shown that there is a unique solution to the problem in our setting). 


Remark 41. This ik is no longer the same as in the previous case, but we 
keep the notation to he consistent with the original artieles as our arguments 
rely heavily on the results in the original articles. 


Definition 42. For each s € (—iM+) (i.e. purely imaginary complex number 
with strictly negative imaginary part), let ((i—^(C, s) he the unique solution 
to the following Riemann-Hilbert problem: 

Find a funetion ik = 'k(C, s) such that 
• 4':C\r^C2x2 is analytic, where 


(148) 


r = Ti = i r2 = i + e* 

r3 = -z + e*^M+, r4 =-i + e®^]R+, r5 = [-i,i]. 


• ik satisfies the jump conditions: for C S 


(149) 4/+(C,s) = 4/_(C,s)Jfc, 

where \k_|_ (^-) denotes the limit of 4/ from the left (right) of the 
contour (the arrows in Figure 1 determine the orientation of the 
curves), and 


(150) 


Ji 


J 3 


J 5 



J 2 


Ja 


(J 'T) 


• In all regions. 


4/(C, S) = {I + 4/i(s)C-' + 'k2(s)r' + 0(|cr')) 

as C ^ 00 . 


( 151 ) 
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• For f3 ^ let 


(152) g = h = - \ (i-e--/^) = ge--^ 

2^sln(7rp) 2zsin(7rp) 


(153) 


Gin = ( J 1) > Gi = GniJt^ Gn = GiJi, 


and 


(154) Hjii = ^ , Hiv = HiiiJ^^, Hi = HivJi^. 

Then define Fi = Fi{(,s) in a neighborhood of i by 


(155) ^{C,s) = FiiC,s)iC-i)^^^G, 

in regions j = I, II, HI, where the branch cut of (( — i)^^^ is along 
i + e~ (0, oo) and arg((( — i) G (—57r/4, 37r/4), and F 2 = F 2 {C, s) in 
a neighborhood of —i by 


(156) ^iC,s) = F2{C,s){C + i)^'^^Hj 

in regions j = I, III, IV, where the branch cut now is along —i + 
e~ (0, 00 ) and arg(C + i) G (—37r/4,57r/4). 

If (3 & define Gm = Hm = I and the other matrices through 
the jump matrices as before, and define in the region j, Fi and F 2 
(with similar branch cuts) by 

(157) 4/(C,s) = Fi(C,s)(C-i)2"^ G,- 

and 

1 log(C + *) j pj. 

Then these functions Fi and F 2 must be analytic functions of C 
in some neighborhoods of hi. 

We thus find 


^Ar(^) = -'^ 11 ( 2 ;) 


(159) +Ti2(z) 


Di{ft,z)Deift,z)\^ Z2 




'^21 ( - log z, —2iNt 


1 JV 

2 2 2 


'I'll - log 2, —2iNt 


.Drift, z)Dfift,z)J ^/i(z) \t 
We note that if we write s = —2iNt and C, = j log 2;, then 
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Figure 1. The jump contour and jump matrices for 'h (a 
modification of Figure 1 in [6]). 


(160) =64*^ = 641^1^ 

so the relevant question is the asymptotics of the hrst column of 

(161) ^'(C,s)e3l*l^^3^ 

where s is purely imaginary with a negative imaginary part and as we con¬ 
sider z = {1 — r)e**, C = * + j log kl- We need to study these asymptotics 
in the different regimes of C, and s (C close to i or far away from i; s small, 
bounded, or large). This can be done by making use of results in [6]. 

Remark 43. For notational simplicity, we will only go over things for non¬ 
integer {3. Integer fl would eomplieate some formulas slightly, but can be 
done in a similar manner. 

Bounded s: Let us first focus on the case where s is bounded, i.e. we 
assume that there exist some constants 0 < c < C < oo (we think of c being 
small and C large) and consider c < Nt < C so that s = —2iNt is bounded. 
Let us now assume that is in a small neighborhood of i so we have 

(162) 4^(C,s)e^'^^ =Fi(C,s)(C-i)="^Gme^"L 

As Gjii is upper triangular with bounded entries and Fi{(, s) is bounded 
in the neighborhood of i (as s is bounded), we see that for in the fixed 
neighborhood of i, 
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(163) 


klc 

^(C,s)e—“"s 


1^0(l)(C-i)^ *) 


where the 0(1) terms are uniform throughout the neighborhood and uniform 
in c < iVt < C. 


If on the other hand Q is outside of the given neighborhood (and as we are 
outside the fixed neighborhood of —i), we have by the asymptotic expansion 
of 'h - (151) - that 


(164) 4'(C,s)e^^^3 = 0(1) 

uniformly in the relevant parameters. 

Let us now note that for G [1 — 7 ,1 + 7 ]e**, Q — i = (log Thus Q 

in a small neighborhood of i means roughly (recall t N-^) iV|l-|z|| < 1. 
We conclude that 


(165) 


0(l)^(iV|l-|z||)2, iV|l-|z||<l 

= { V 

^O(l)^^, N\l-\z\\>l 


VhD ’ 

uniformly in 2 : and Nt G (c, O). 

Again the relevant question is to analyze the first column of 


(166) 






Let us start with the case where C is close to i (so |1 — I^H < W ^). Let 
us make use of the representation 'L((^, s) = s)(( — We have 

(denoting by T’{(^, s) the derivative of Fi((, s) with respect to () 


d 


(167) 




VM- 




+ ^T((C,s)(C-i)2"3Gme^^"3 

+ A^^Fi(C, s)^3(C - 

zt Q — i 2 


For the first term, we note that 


(168) ft{z) = (1 - 2 e**)t(l - 2 e-**)t(l - 2 -'e**)t(l - ^-^6"**)^ 


so that 
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( 169 ) 


fliz) 


/3 


ft{z) 

We then point ont that 


z — e‘' 


+ 


/3 


z — e 


wr + 0(i). 


(170) 

(Fi(C, s)(C - = (Fi(C, s)c73(C - , 

so we can combine the first and third term of (167) into 


(171) 0(1) 






/3 1 


‘1 z — e* 


/3 ^ 1 1 


z — e 


2 zt C — i 


Writing l/((^ — ^) = t/i\og{z/e'"^)) we see that this can be written as 


(172) 


0 ( 1 ) 




(C 



= 0 ( 1 ) 


1 


iV(iV|l - 


+ 0(l)t-i + 


u -1 


1 


\z\ log 1^1 


^11) 


i. 

2 


with 0(1) being uniform in z and Nt G (c, C). The remaining two terms in 
(167) can be estimated to be of the same order so we find that for Ai|l —12:|| < 
1 , 


(173) T)v(^) = 0(l)^^iV(iV|l - |z||)t. 

For |1 — l^ll > N~^, we note that it is proven in [6] that 

(174) ^^{C,s) = AiC,s)^{C,s), 

where for bounded s, A is (uniformly) 0(1) outside the fixed neighborhoods 
of ±i. Thus 


(175) 

and 




= 0 ( 1 ), 




iT(C, s)e = 


1 - Ul 




0 ( 1 ) 


(176) 


= 0(1)W 




Large s: Consider now the situation where Nt > C. For this, we make use 
of the s —>■ —ioo asymptotics of 'k(C, s) which have been studied in section 5 
of [6]. For us, the central result of this section is the following: in the region 
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we are interested in (i.e. for z = |z|e** with \z\ < 1 and small perturbations 
of the angle), 


[R[C,s), IC - > 1 

where R{(^,s) = 0(1) and ■^R{C,s) = 0(1) uniformly in ^ off of a certain 
jump contour which we don’t need to worry about and uniformly in the 
relevant domain of s. 

Remark 44. Here in the condition |C ~ < 1? the upper bound should 

be replaced by some small enough 0(1) number, but this only changes our 
results by a uniform 0(1) factor, so for notational simplicity, we stick to 
this condition. 

Here the function Pi is defined in the following manner: 

(178) Pi(C, s) = (^^(C - i)^ 

where in the region relevant to us, 



(179) M(A) = 

where L is a function which is analytic in a neighborhood of zero and depends 
only on (3. The branch of the power is chosen so that argA G (0, 27 r). G3 
is an upper triangular matrix with ones on the diagonal. Another fact we 
need about M is that as outside of a neighborhood of the origin. 


(180) M(A) = (1 + 0(A-i))e-5^'"T 

Let us consider first the situation A^|l — |2;|| < 1. This translates into (up 
to a uniform 0(1} factor) |s||(C — z)| < 1 so as |s| is large we have 


T(C,s)e^^'"3 = 0(l)(|s|(C - 

(181) = 0(1)(A^|1 - 

_ (0{l)(N\l-\z\\)2e^^ *\ 

VC>(l)(iV|l - |z||)2e^^ *)’ 

We then note that 

(182) = = 0(1) 

for the f we are considering now. We conclude that for A^|l — \z\\ < 1, 


(183) 


^n{z) = 0(1)^^(N\1-\z\\)^ 

VftG) 


uniformly in the relevant variables. 
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Consider next the situation where A^|l — \z\\ > 1, but |C ~ < 1 (these 

conditions are roughly equivalent to N~^ < |1 — |^;|| < t). In this region we 
use (180) for estimating M. We find 

(184) ^(C, s)e^^'"3 = 0(l)e-^(^-*)e-x/3^3g^C^3 ^ 

and conclude that 


(185) <^n{z) = 0{1)^^ 

uniformly in the relevant parameters. 

Consider finally the case where |C ~ > 1 (or |1 — lz|| > t). Here we have 

I s I - 

also 'I'(C, s)e^'’°'^ = 0(1) so again 


(186) <^n{z) = 0{1)^^ 

uniformly in the relevant parameters. 


For the derivative, we note that in the \\z\ — 1| < N~^ regime, the ar¬ 
gument of the bounded s case goes over in an analogous manner (so dif¬ 
ferentiating picks up a factor of \s\/t = 2N) and we have in this regime 
(uniformly) 


(187) 4>)v(^) = 0(l)^^iV(iV|l - \z\\)l 

In the regime N~^ < |1 — I^H < t, ^('I'(C, s)e^‘’^®) = 0(1) so one finds 
with similar arguments 

(188) =0{l)^^\l-\z\\-\ 

Finally for |1 —12;|| > t we have again with similar arguments for estimating 
ftiz) and its derivatives 


(189) ^'^{z) = 0{l)^^t-\ 

Again all of these estimates were uniform. 

Small s: Finally we consider the case where Nt < c. For this case, it is 
proven in [6] (section 6) that if we write A = |s|(C — i)/2 = A^log \z\, then 
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( 190 ) 




H{X,s)Po{X,s), iV|l-|z||<l 
H{X,s)M{X), iV|l-| 2 ||>l 


where M is as in (178) except that (3/2 is replaced by j3. In the region 
relevant to us, Pq can be written as 


(191) Po{X,s) = L{X) 

with a suitable convention for the branch cuts (these aren’t important for 
us as we are estimating only magnitude). The function L is analytic in the 
relevant neighborhood of the origin (in fact it’s essentially the same function 
as before). The function J can be defined as 


(192) J(A;s) = - 

Js ? - A 

Also H and its derivative are 0(1) uniformly off of a certain jump contour 
(whose details are not relevant for us). 

Remark 45. Again the neighborhood of the origin is not necessarily exactly 
A^|l — |z|| < 1, but up to a 0(1) factor, we can pretend it is. 

Consider then the A^|l — |z|| < 1 situation. Arguing as before (e.g. the 
(1,2) entries of the upper triangular matrices don’t affect the first column 
of T), we find 


(193) 


'I'(C, s)e _ 


/o(l)(A^log |z|) 2 (A^log |z| — s) 2 
I 0(l)(A'log l^l) 2 (A^log |2:| — s)^ 




e 4 




e 4 


We hrst note that A^ll — |2:|| < 1 is roughly equivalent to |s||C — < 1 so 

the last two terms don’t affect anything. We then note that 


(194) (Wlog l^l — s )2 = 0(l)Ai2 (log l^l — 2it) 2 
so we see that 

(195) 4>7v(. 2) = - A^^|l — |z|| 2 max(|l — \z\\i,t^). 

For 77|1 — l^ll > 1 we recall (180) and hnd (uniformly) 

4'(C,s)e^‘’°'’^ = C)(i)e“i^^°§l^l'^3g^0+jiogb|)cr3 

(196) =0(1). 

Thus for A^|l — l^ll > 1 

(197) ^n{z) = 0{1)-^ 

Vftiz) 


uniformly. 
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Let US now consider the derivative. Returning to the notation A = |s|(C ~ 
z)/2 = A^log \z\, we recall that we can write (for A^|l — \z\\ < 1 or roughly 
equivalently |A| < 1) 


(198) 


dl'(C, s)e _ 


f0{1)X2 {X - S)2 

\C>(1)A2(A-s)2 


where the 0(1) also have uniformly 0(1) A-derivatives. We then note that 


(199) 

We thus have 


d Id 1 |s| d N d 

dz tz dQ tz 2 dX z dX 


d 


1 


wCo-3 


( 200 ) 

1 


^(C,s)e 4 


'0(l)At(A - s)t (-l|g + f f (A-i + (A - s)-i) 


V \C)(1)A2 (A - s) 2 (Y 2 ^ “'S) ^) 



N fo(l)X 2 (X- s )2 
VM^) \0(i)X2(X-s)2 
Recalling (169), we have 


sfi -a"- 


/3 1 




2 z — e* 


2 z — e 


—it 


+ 


I3N 


1 


2z Nlog \ z\ 


( 201 ) 

( I 1 ^ 

2 \1 — l^l \z\ log \z\) 

= 0 ( 1 ) 





/3N 1 

2z A^log |2;| + 2iA^t 


+ 0 ( 1 ) 


1 1 1 \ 

26** — 1 ze** log(2e**) J 


+ 0 ( 1 ) 


uniformly in the relevant parameters (here we used the fact that ((1 —r(;)“^ + 
l/(t(;logic) is bounded for w close enough to 1. We conclude that the 
relevant asymptotics of ^']si(z) for A^|l — |z|| < 1 are 


cl>)v(^) = 0(l)iV^^Ai(A-s)f 

Vfto) 

(202) = 0(1) ^ 1V^+^|1 - \z \\2 (\og\z\ + 2 it )2 

Vftiz) 

= 0(1)— \z \\2 max(|l — jzll + tf) 
^ft(z) 


uniformly. 
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For A^ll-^l — 1| > 1 we had 


^(C,s)e^^'"3 = 0 ( 1 ) 


(203) 

uniformly as is its A-derivative so we find 


yfti 


VM^) V 

(204) = 0{1)N 

We conclude that for A^|l — jzlj > 1, 


(o(l)N + Oil)jl^ 




(205) 


\/ft(z. 


□ 


We will also need a result on the asymptotics of (pwiz) and in the 

case where Nt < c, but is not on the contour, but in the sector between 
(1 — e, 1 + e)e^**. These asymptotics also follow directly from the analysis 
in [6]. More precisely, 

Lemma 46. For Nt < c, where c is as in the previous theorem, and z = re*®, 
where r G (1 — 7 , 1 + 7 ) and 9 G (—t, t), we have: for A^|l — |z|| < 1 


<hAr(z) = 0(1)— ^ - max(|l — \z\\'^ ,{t — 0)2)max(|l — | 2 ;| 12 , (t + 0) 2 ) 

VM^) 

(206) 

= 0(l)iV^ 

and 


^' n ( z ) = 0(1)-^^-max(|l — 1^11 2 , (t — 0 ) 2 ) max(|l — |z|| 2 , (t + 0) 2 ) 

\/ft(z) 

(207) = 0(l)W^+i 
For N\\z\ — Ij > 1, 

VftiN 

(208) = 0(1)11 - |z|r^ 

and 

4>)v(z) = 0(l)iV^^ 

(209) = 0(1)A^|1 - |z|r^. 
uniformly in z and t < cN~^. 
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Proof. As in the previous proof, let us write C = log z/t and s = —2iNt. We 
also consider for simplicity the \z\ < 1 case. Again, the relevant quantity is 
the first column of 


( 210 ) 


and the ^;-derivative of this. Let us consider first the case where z is close 
to or A^|l — \z\\ < 1. We hnd as before (with A = |s|(C — f)/2) that 


( 211 ) 


|s| - 

'L(C, s)e 4 


/c>(1)A2(A-s)2 *\ 

\0(l)A2(A-s)t *J ' 


We note that in our case. 


(212) \ = '-^{C - i) = Nlogz - iNt 

and 


(213) A — s = A^log 2: + 

We conclude that 


A2 (A — s) 2 = C>(1)A^^ max(|l — \z \\2 — 6 ) 2 ) max(|l — \z \\2 ^ {t + 6 ) 2 ). 

Noting that in our case 

^^ = 0(l)|2 -e*r^k-e-*r = 

(214) =0(1) min(|l — |2;||“2 ^ (f — 0)“2) min(|l — |2;||“2 ^ (f + 0)“2) 

so our claim about the asymptotics of $7v is immediate. 

For iV|l — l^ll > 1, we have 


(215) ^j,{z) = 0{l)—^ 

For the derivative in the A^|l — |2;|| < 1 case, we saw that 
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d 


1 


|s| - 

V^°'3 


V ftiz) 

(216) 

1 


^(C,s)e 4 


^0(l)At(A - s)t (-i|g + f |(A-i + (A - .)-!)' 


V \C)( 1 )A 2 (A - s) 2 (^-1A^ + :^|(A ^ + (A - s) 


X 2 ' 


^ iv /o(1)A2(A-s)2 

v//t(2:) yo(l)At(A - s)2 


As earlier, one sees that the second term is dominant and differentiating 
picks up a factor of N. 


(217) 


Consider now the A^|| 2 :| — 1| > 1 situation. In this case we have as before 

1 


^)v(^) = 0(l)iV 


Vm- 


□ 


5. Asymptotics of the Fredholm determinant 

Our approach for analyzing the asymptotics of the Fredholm determinant 
is based on the following estimate (see [30] Theorem 6.5 - set B = 0 and 
A = K). 

Theorem 47 (Simon). For a trace class operator K, 

(218) I det(/ + K)e-^^^ - 1| < || ATHae^^dl/flh+iF 

for a suitable positive constant V 2 , where the determinant is a Fredholm 
determinant and \\K \\2 denotes the Hilbert-Schmidt norm of K. 

Thus to estimate det(/ + K), we must be able to estimate ||Ar ||2 and trAT 
(in fact we’ll conjugate AT by a certain multiplication operator - the trace 
and the Fredholm determinant remain invariant, but the HS-norm does not). 
Let us do this in the three different cases. 

5.1. No singularities. We now make use of the asymptotics of the orthog¬ 
onal polynomials from Theorem 34 to estimate the asymptotics of K in the 
case where the symbol has no singularities. This was essentially done already 
in [3], but for completeness (and to make sure of the relevant uniformity) 
we give an argument here. 

Lemma 48 . Consider a Laurent polynomial V of the form 

(219) F(z) = -doj-z"^). 

j =0 

Also fix some R > 0 and write Aji = {maxj jajj < A}. Then there exists 
a a > 0 such that as N,M ^ 00 (such that M > N), 
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(220) det(/ + K)= 
uniformly on Ar. 

Proof. As noted in Remark 32, as the symbol is analytic in C \ {0}, the 
integrals across the rays will cancel and we only need to worry about the 
behavior of the kernel on the circles (1 ± e) x T. 

We also note that as the Fredholm determinant is invariant under conjuga¬ 
tions (i.e. for invertible X, det(X(/-|-i^')X“^) = det(/-|-iF)), we can replace 
K by where : L'^([JjTj,dzl2Triz) L‘^{UjTj,dz/27riz), 

{AI]\ih){z) = z~^^'^h{z). This effectively symmetrizes the operator K. The 
kernel of this operator (which still is trace class) is {z/w)^^'^K{z,w). 

Consider first the case \z\ = 1 — e, |u;| = l-|-e (or vice versa). We have 
by Theorem 34 (uniformly) 

N _ N _ 

(7) " - ii;) " 

_ W 

Z 

(221) =0(l)(l + e)T(l-e)-f. 

Thus 


(-)" Kcd{z,w) = 
\w/ 


(222) K{z,w) = 0(l)(l + e)“'^(l-e)'^(l + e)^(l-e)"f = 0(e“"(^"^)) 
for a suitable a > 0. 

Consider then the case \z\ = |tc| = 1 ± e. In this case, we write 

N _ N _ 

(7) " - (7) " 


(-)^ Kcd{z,w) = 
\w/ 


1 - 


(223) 


w^(t)p^{w ^)-z^(t)j^{z ^) 


= z^ 2 re 2 (pj^{z) 


NZ 


z — w 


, 1 +^ -N- ]^ f)N{z)-(l)N{w) 

+ Z^2w ^(Pn[z )- 


Z — W 


We then note that for a nice enough function h. 


(224) 

for some point ^ with |^| 


h{z) - h{w) ^ 
z — w 

\z\{= Iml) so we see that 


Kcd{z,w) =0{l)N(j)N{z)(t)j^{C, ^)PO{l)(t)N{z)(l)M{C ^) 

(225) +O{l)^^{z-^)cf'^{0 

for some points with |(^| = |^| = \z\. In each of the terms above, at least 
one of the two points is on the circle (1 — e)T giving exponential smallness. 
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This exponential smallness cancels the factor of N, so we see that all of these 
terms can be bounded (uniformly) by 0(1) max(|z|'^, We conclude 

that for some a > 0. 

(226) 

(' K{z, w) = 0(1) min(|z|^, l^l"^) max(|2|^, l^l"^) = 

\wJ 

uniformly in z and w as well as on A^. 

These estimates for the kernel imply that 

(227) tri^, IIMtv1 12 = 0(e-“(^-^)) 

uniformly on A^. Making use of Theorem 47, this gives the desired result. 

□ 

5.2. Distance between singularities bounded away from zero. Mak¬ 
ing use of Theorem 36, we try to mimic the proof of the non-singular case to 
analyze the asymptotics of det (H-iiT) in the case when the distance between 
the singularities is bounded from below. 

Lemma 49. Let 17 : C \ {0} —>• C 6e a Laurent polynomial of the form 

(228) V{z) = +a]z-^), 

j=0 

{wi, ...,Wk) be distinct points on the unit circle such that 0 < arg(t(;i) < 
arg{w 2 ) < ... < aicg{wk) < 27r, and assume that there exists some fixed 6 > 0 
such that \wi — Wj\ >6 for i 7^ j. For \z\ = 1, let 

k 

(229) f{z) = e^^^^ll\z-WjfL 

1=1 

where (dj > 0, and for any fixed i7 > 0, Ar = {maxj \aj\ < ii}. Assume 
further that there exists a q € {0, 1) such that for large enough N, N/M < q. 

In this case, 


(230) det(/ + i^) = 1 + 0(1) — 

uniformly on Ar and uniformly on {miuj^j \ wi — Wj\ > 6}. 

Proof. As in the non-singular case, let us conjugate K by the multiplication 
operator Mjv- We also start with the case \z\ < 1 — 5/6 and |t(;| > 1 + 5/6 (or 
vice versa). Making use of Theorem 36, we have (uniformly in everything 
relevant) 


(231) 



Kcd{z,w) 


0(1) maxdzl 2 , |2:| 2 ) x max(|r(;| 2 dtc| 2) 


A DISCRETE LOG GAS, DISCRETE FH-SINGULARITIES, AND GMC 


45 


and for some a > 0 


(232) (-) " K{z,w) = 

\w/ 

The rest of our analysis we split into two cases: \z — w\ < N~^/2 and 
\z — w\ > N~^/2 (we call these the short and long range regimes). The 
point of this is that for \z — w\ < N~^/2, we replace difference quotients 
by derivatives and for any point C, with jz — C| = = 0{\)z^. 

For \z — w\ > N~^/2, we simply replace the denominator 1 — z/w in the 
CD-kernel by N~^. 

The short range regime. As in the non-singular case, we can write 

/ z\- - 

(-) ^ Kcd{z,w)=0{1)N4>n{z)^n{C^) + 0{1)cI>n{z)4>n{C^) 

\w/ 

(233) +0{l)4>j^{z ^)</>)v(C) 

for some points C between z and w on the contour (here we made use of 
the fact that for \z — wl = 0{N~^), |C|^, = 0{l)\z\^ uniformly 

in everything). 

If \z\, |t(;| < 1 — 5/6 or \z\, |r(;| > 1 + 5/6, then by Theorem 36, we see that 
reasoning similarly to the non-singular case, we obtain the uniform bound 

^ Kcd{z, w) = 0(1) maxd^l^, I^T^) 

\wJ 

H. —H. H. —H. 

(234) = 0(1) max(|z| 2 ^2:1 2 ) x max(|t(;| 2 , |t(;| 2 ). 

This then implies that for some a > 0, we have uniformly in everything 


(235) K{z,w) =0(e-^^^-^'>). 

As the distance between the singularities is bounded away from zero, this 
only leaves the case z,rt; G [1 — 5/6,1 + 5/6] x Wj for some j (the case 
where one of the points is near the end of the interval and one right next 
to the interval can be treated as the case where neither is in the interval by 
retuning the parameters). Here we need to note that as jzj —1, we don’t get 
exponential smallness from v(z). In fact, in the worst case scenario, where 
Wj G Vm (i-e. = 1), v{wj) = 00 . We assume the worst case scenario is 

the generic one and we use the following estimate for v: 

0(1) min(|2:|^, |2:|“-^), |1 — |2;|| > 

Due to this change of behavior at a scale of M~^, we need to split our 
treatment further into several subcases. We also will need to make use of 
the fact that f{z) = 0(1)|1 — \z\\^^ is present in the kernel. 
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|1 — \z\\, |1 — |u;|| > N In this case, using Theorem 36, we have (note 
that in this case, we have for example |1 — |C||~^ = |1 — |u)| + |r(;| — |C1I~^ = 
0(1)|1 - |u;||-^) 

(237) (-)^ Kcd{z,w) =0(1)A^|1 - |z|r^|l - |u;|r^ 

\w/ 

so that 

(238) ( — ) ^ K(z,w) = 0(1)N min(| 2 ;|^, bl”^) min(|r(;| ItcT^). 

\w/ 

Again by retuning the parameters, we see that essentially the only other 
possibility for |z — u;I < A^“^/2 is that |1 — |z||, |1 — |t(;|| < In this case, 

Theorem 36 implies that 

(239) (-)~ Kcd{z,w) =0{l)N^+f^i 
\wJ 

uniformly. One then needs to take into account the different behavior of 
v{z) for |I — |z|| < M~^ and |I — |z|| > M~^. We have the following cases: 

M~^ < |1 — l^ll, |1 — |t(;|| < N~^\ Here we have 


/ Z \~2 / X /'/-i/x /iiMti — Mx /I iMt \ — M.\ 

— K{z^w) = O(I) mini\z\ 2 Az\ 2 ) min(ki? 2 1^; 2 ) 

\w/ 

(240) X |1 — |z|k|l — 

|1 — l^ll < M~^ and M~^ < |1 — |r(;|| < N~^ (or vice versa): In this case, 
we have 

(241) 

^ ^ K{z,w) = 0(1)11 — \z\\'^~^\\ — \ w\\^ min(|rt;|~, \ w\~~) 

and the other case by interchanging the roles of z and w. 

|1 — |z||, |1 — |t(;|| < M~^: Finally we have 


(242) 


z\f N^+0j Pj 1 Pj 1 

-)' K{z,w) = o(i)^^ii - \z\\^-m - kii^-2 


w 


M 


uniformly. 

The long range regime: In this case we estimate 


N 

Z\^ 


w 


( 243 ) 


Kcd{z,w) = 0{l)N 

+ 0{l)N 


N 

2 \ T 




W 

N 

U)\ T 


(j)Niw)(j)^{z 


4>Niz)4>N{w ) 


Again we consider the different parts of the contour z and w can lie on. 
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zj, Iwl < 1 — 5/6 or \z\, |ti;| > 1 + 5/6: Making use of Theorem 36 we have 


(244) — Kcd{z,w) = Olljuiaiicilzl 2 \w\ 2 Az\ 2 \w\ 2 ) 

\w/ 

and we have as before, for some fixed a > 0 


(245) f-) ' K{z,w) = C)(e-“(^-^)) 

\w/ 

uniformly. 

Let us now focus on the situation where z G [1 — 5/6,1 + 5/6] x rcj and 
w G [1 — 5/6,1 + 5/6] X Wi for some i and j (note that possibly i = j). 

|1 — jz]], ]1 — Ircj] > N~^: In this case our worst case scenario is where z 
and w are on the opposite sides of the unit circle and we have 

(246) 

/Z\~2 /X / X /, M-N M-N M-N M-N 

— Kiz^w) = 0{\)N m.m.{\z\ 2 \z\ 2 )min(kt; 2 \w\ 2 ). 

\w/ 

|1 — jzj] > N~^ and M~^ < jl — jrc]] < (or vice versa): In this case, 
= 0 (1) so we see that 


{—\^K(z,w)=0{l)Nva.m.{\z\ 2 ^\z 
\w/ 


(247) 


li¬ 


re 


' t „ - ' t 

2 N 2 


M-N /, I M I 

2 ) min(|t(;| ^ ^\w 


M 

2 


) 


The opposite case is obtained by interchanging z and w. 

|I — | 2 ;|] > N~^ and jl — jte]] < M~^ (or vice versa): In this case, we need 
to take into account the behavior of v near the unit circle. We find 


) — ]\T 

^ K(z,w) = 0(1)—;= min(jz|“ 

Vm 


,M-N M-N . 

U 2 


(248) 


I I 11^_L AT — 

X 1 — ke 2 2N 2 . 


The opposite case is obtained by interchanging z and w. 


M ^ < ]1 — jz]], jl — jre]] < With similar arguments. 


- 1 . 


( 249 ) 


N 

2: \ T / N / N /, I M I— Mn |M 

—j K(z,w) = 0(1)N miu.(\z\ 2 ^\z\ 2 ) min(jrel 2 , |u;| 2 

ft ft ft ft 

X 1 - ke 1 - \z\\-N-. 


M ^ < jl — \z\\ < N ^ and jl — |r(;|| < M ^ (or vice versa): Here wehnd 
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(250) 


N TKT 

— K{z,w) = 0(l)^=m.ui{\z\‘2 ,\z\ 2 ) 
wJ 'Vm 

Pi 1 Pi 

X 11 - \w\\^-2N^\1 - Ull^iV^. 


The opposite case is obtained by interchanging z and w. 
|1 — 1^11, |1 — Iwll < M~^: Here we have 


(251) 


N 


w 


N 




1 


- K{z,w) = 0(1) —11 - \z\\^~^\l - |u;||^- 2 iV^. 


M 


We are now in a position to estimate the trace and the HS-norm. Let us 
start with the trace. 

Estimating the trace: As K is trace class, its trace is well defined and 
given by K{z, z)a{dz). Thus we only need the short range asymp- 

totics {z = w). We then decompose our integration contour into the relevant 
parts: 


/-/ "/ 

JTi Jrwnlli-kll><5/6l Jr 


(252) 


n{|i-L||><5/6} ./rjn{Af-i<|i-b||<5/6} 


+ 


By (235),(238),(240), and (242) we have then for some fixed a > 0 


trA: = 0(e-"(^-^)) +V / 0(1)A^(1+ r)-^dr 

+ V / 0{l)Nf^^+\l+r)-^rf^^dr 

+ E 


/ 0(1)—77 —^dr 

lo M 


(253) 


N 


N 


= Ole-X"-")) + 0(1)- + 0(1)- + 0(1) ( - 


M 


-) 

M J 


i+^j 


= 0 ( 1 ) 


N 

M 


uniformly on Aji. 


Estimating the square of the HS-norm: To apply Theorem 47, we 
wish to also estimate the HS-norm of . By definition, this is given 

by 

(254) \\MmKM]^% = EE// (T((i 2 :)cr(dr(;) I — I \K(z,w)\^. 


j I 
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We decompose Fj x F; into the short and long range regimes and then 
split these into parts where we have estimated the different asymptotics of 
the kernel. Let us consider first the short range regime. We split Fj x F; 
in the short range regime into the cases where i) |1 — |^;1|,|1 — |tc|| > (f/6, 
ii) |1 — |z|l,|l — \w\\ G {N~^,5/6), hi) |1 — |^;||, |1 — |tc|| G 
iv) |1 — l^ll < M~^ and |1 — |t(;|| G {M~^, N~^) (or vice versa), and v) 
|1 — \z\\, |1 — |rc|| < M~^ (recall that the remaining cases can be included 
into these by retuning the parameters). Making use of (235), (238), (240), 
(241), and (242) we see that the contribution of the short range case to the 
square of the HS norm is 


+ 0{1)N^ \^j (1 + r)-^drj 

+ 0(1) ^ (1 + r)~^r^^dr^ 

^dr / (l+r)-^r^idr 


Ar2+2/3j rM 

(255) +0(l)Y,^^j 


+ 0(1)E 


^2+2/3j / /-M-l 
M2 


I 


^Pj-^dr 


= 0 ( 1 ) 


N 

M 


For the long range regime, we again split it into the relevant subregions: 
i) |1 — \z\\, |1 — \w\\ > h/6, ii) |1 — | 2 :||, |1 — |rc|| G (N~^,S/6), hi) |1 — \z\\ G 
{N~^,6/6) and |1 — |ic|| G (M~^, N~^) (and vice versa), iv) |1 — \z\\ G 
(W“^,h/6) and |1 — |r/;|| < M~^ (and vice versa), v) |1 — |2:||,|1 — |rc|| G 
(AI~^, N~^), vi) |1 — \z\\ G (M~^, N~^) and |1 — |rc|| < M~^ (and vice 
versa), and vh) |1 — | 2 ;||, |1 — |rc|| < M~^. Making use of (232), (245), (246), 
(247), (248), (249), (250), and (251) we see that the contribution of the long 
range regime to the square of the HS-norm is (for some fixed a > 0) 


0(e""(^"^))+ 0(1)A^2 ' (i + r)-^^-^^dr^ 

rS/6 

+ 0(1) / (l + r)-(^-^)dr / (l + r)-^r^Mr 

, iw-i 7 m-i 


3 

AT^+di /■<5/6 /-M-i 

+ 0(1)--—Y (l+r)-(^-^)dr / rf^^-^dr 
A'l ^ Jn-^ Jo 


M —JM 
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(256) + 0{1)'S^ / {l + r)~^r^^dr / (l + r)“^r^*dr 

^ Jm-1 Jm-1 


?Y2+/5i+/5_j' z*-^ ^ z*-^ ^ 

+ 0 ( 1 )^^-—- / {1 + r)~^r^^dr / r^"~^dr 

^ Jm-1 Jo 

^iJ 

7Y2+/3i+/3j 

+ 0(1)^-- / r^^- 


M2 




/ r^^~^dr / r^'~^dr 

Jo Jo 




where we made use of N/M < q < 1 for large enough n. 

Combining these, we see that = 0{1)N/M so that by 

Theorem 47, 


(257) det(/ + i7) = 1 + 0(1) — 
uniformly on Ar and {minj^j \ wi — Wj\ > <5}. 

□ 

5.3. Merging singularities. The basic idea here is the same as before. We 
estimate the kernel with the help of Theorem 39. For simplicity, we use the 
symbol \z — e^^\^\\z — . The estimates are the same in the case obtained 

by rotating this symbol. 

Lemma 50. For \z\ = 1, let ft{z) = \z — Assuming that 

there exists a q G {0, 1) such that for large enough N, N/M < q, we have 

(258) det(/ + i7) = 1 + 0(1)^, 
where 0(1) is uniform in 0 < t < to- 

Proof. Let us begin by considering what is different in the case of merging 
singularities. First of all, the asymptotics of the orthogonal polynomials are 
different. Next we note that if t = 0(1)we might have the issue that in 
the short range regime, we need to consider the case where 2 ; € (1~C) l+e)e** 
and u; e (1 — e, 1 + e)e“** (or vice versa). In this case, we can’t use the same 
estimate for the CD-kernel we used before - the points C and ^ might not 
be on the contour any more (or then they might not be close to z and w). 
For this, we will use Lemma 46. Another issue is that the asymptotics of ft 
near are different. Instead of f{z) = 0(1)|1 — \z\\^, we have 

(259) 

/^(z)=0(l)|l-|^||^||z|-l + (l-e-2**)|^ = 0(l)|l-|^||^nax(|l-|z||^^^). 

Apart from these differences, the reasoning is similar to the previous case 
so we will not be as detailed as we were there. We will again consider the 
conjugated operator Mat ATLet us begin with considering the long 
range regime. 
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Long range regime: As in the case of non-merging singularities, we 
estimate the CD-kernel as 

JV JV 

(w) " ^cd{z,w) = 0{N) ^ (/)jv(z)'^n(w~^) 

(260) +0(N)(-)^^jy(z-^)Mw). 

\w/ 

As the polynomials have the same asymptotics as in the case of the non¬ 
merging singularities for |1 — |z|| > 7, we can focus on the case where at 
least one of the z, w is at distance less than 7 from the unit circle. We then 
recall that K{z, w) comes with a factor of -\/ ft{z)-\/ft{w) so what is relevant 
is the asymptotics of ^^fti^cpNiz)- Comparing Theorem 36 and 39, we see 
that these are the same in the two cases (the only case that looks different 
is the A^|l — l^ll < 1 case for Nt < c, but here the difference is the term 
max((A^|l — jzll) 2, (A^t) 2) < max(l,c2)). Thus in the contribution of the 
long range regime to the square of the HS-norm of is the same 

as before (of order (N/M)'^). 

Short range regime: We again estimate 

(-) ^ Kcd{z,w)=0{1)N4>n{z)^n{C^) + 0{1)cI>n{z)4>n{C^) 

\w/ 

(261) +0{l)4>j^{z ^)</>Ar(C) 

for some points ( with |z — C| = 0{N~^) and |z — ^| = 0{N~^). If z and 
w are on the different rays, we take ^ and () such that they are on the line 
going through z and w (this can only happen if t = 0(1)A^“^ which up to 
a 0(1) term is equivalent to the Nt < c case). Otherwise, we take them on 
the contour and between the points. 

For |1 — l^ll > 7 and |1 — |rc|| > 7 we can reason as before to get expo¬ 
nential smallness. Let us then simply argue as in the case of non-merging 
singularities to find the asymptotics of the kernel for different values of z 
and w. 

|1 — \z\\, |1 — lirll G (iV“^,7].- We have now e.g. |1 — \z\\ = 0(1)|1 — |C|| 
implying also ftiz) = 0{l)ft{C)- Let us first consider the case where z, w 
are on the same ray (as the one where they are on different rays is only 
relevant for Nt < c). Using Theorem 39, we find 

(262) (^) ^ Kcd{z, w) = 0{1)N 
and 



N 

T M 

K{z^w) = 0{1)N min(|2:| 2 




M 

2 


min(|i(;| 2 


w 


M. 

2 


) 


(263) 


w 
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uniformly in z,w and 0 < t < Iq. Again up to a uniform 0(1) term, the 
only remaining situation is when |1 — \z\\, |1 — |u)|| < N~^ and in this case 
Theorem 39 implies that 


(264) 


N 

2\ T 


Kcd(z,w) =0(l)A^^+^(min(iV,t-i))^. 


We thus have (uniformly) 


N 

z \T 


K(z^ w) = 0(l)-\/ v(z)-\/v(w)\l — |z|| 2 |l — |t (;||2 


(265) 


X max(|l — l^ll 2 ,1 2 ) max(|l — |rc|| 2 ,1 2 min(A^^, t ^). 


Consider hnally the case where z and w are not on the same ray (in the 
case where Nt < c). In this case, we estimate (using Lemma 46) 


(266) 


N 

2 \ T 


W 


Kcd(ZiW) = 0(1)N m.m.{\l— \z\\ A^^) min(|l — |t(;|| 


Estimating the trace: For the trace we only need the short range 
estimate. Using (263) and (265), we have for some a > 0 

trA: = 0(e““(^“^)) + 0(l)iV P {l + r)-^dr 

Jn-^ 

(267) +0(1) f |u( 2 ;)||l — l^ll^ max(|l — |z||^, 

The second term is (uniformly in 0 < t < to) of order N/M. So it remains 
to estimate the last term. If t > N it is 

(268) 0(1)N^+^ [ \v(z)\\l-\z\fdz 

which is of the kind of integrals we’ve already estimated and found to be of 
order N/M. If t < N~^, we the term is 


0 (l)^2/3+i /■ \v{z)\\l-\zf^dz 


(269) 


+ 


0{1)N^P+HI3 f \v{z)\\l-\z\fdz 


Here the first term can estimated upwards to 


(270) 0(1)1V2/3 +i /“ \v(z)\\l-\z\f^dz 

./o<|l-L||<Ar-i 
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which again is something we have estimated to be of order N/M. The second 
one can be estimated upwards (recall Nt < 1) to 

(271) 0{1)N^+^ [ \v{z)\\l-\z\\/^dz 

J\l-\z\\<N--^ 

SO we conclude that the trace is of order N/M (uniformly in t). 

Estimating the square of the HS-norm: As noted, the long range 
estimate is the same and of order [N/M]'^. For the short range one, we note 
that if the points are on the same ray, our estimate for the kernel just factors 
into two copies of what we estimated in the case of the trace so from this we 
get a term which is uniformly if order (N/M)'^. We are left with an estimate 
in the case where the points are on different rays (in the case Nt < c). Again 
we have factorization and the remaining estimate is calculating 


(272) 


N^ 


'|i-bll<7 


\ft{z)\\viz)\mm{\l-\z\\-^^,N^f^)dz 


Such integrals were estimated for the trace (and using those estimates we 
find this to be of order (N/M)'^ as well). We conclude that 


(273) ||Mw™-1||2 = 0(1)(^) 

uniformly in t. 

Thus making use of Theorem 47 we have our claim. 

□ 


6. Proof of main results 

Our estimates so far essentially state that if N/M —>■ 0, the discrete case 
does not differ from the CUE. Proving our main results is then essentially 
as in [32]. For completeness (and due to the fact that there are small differ¬ 
ences) we provide proofs here. We will first prove some of the propositions 
presented when describing the structure of our proof and finally combine 
these into a proof of Theorem 9. 

6.1. Proof of Proposition 14 and Proposition 13. 

Proof of Proposition 14- While some variant of the result certainly follows 
from more general results (see e.g. [5]), we choose to give a proof here for 
the sake of completeness. Recall the notation 

N 

(274) Z,=Y.zl. 

k=l 

We also write 


Xj = Re(Zj) and Yj = lm{Zj). 


(275) 
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We wish to prove that for any fixed I, converge in law to 

(Vj72X,-,Vi721S)i=i’ {Xi,Y,,...,Xi,Yi) are i.i.d. standard Gaus- 

sians. Let us consider the moment generating function of let 

ti,Si gU. and consider 

M{ti,ti, si,si) = E 

(276) = ® ( n M(L-“j)4+(h+*s3K'’) 

\A:=1 

If we then define 



(277) V{z) = '^^ {{tj - isj)z^ + {tj + isj)z 

i=i 

we have by Proposition 27 

(278) M{ti, ...,si) = T]\f_i{e^). 

Then by Proposition 30 (and Remark 32) as well as Lemma 48, 


(279) M{ti,...,ti,su...,si) =TN-i{e^) =TN-i{e^) x {1 + 

for some fixed a > 0. Now by the Strong Szego theorem (Theorem 25) 


(280) 


lim T/v-i(e^) = ^ ^ 4^ 

N^OO 


= E 



which finishes the proof. 


□ 


The proof of Proposition 13 is then more or less immediate, and while a 
similar one is presented in [32], we give one for the convenience of the reader. 

Proof of Proposition 13. Noting that 


(281) E(F)^^^_^(e*®)) = (^n f (cos(j0)x,+sin(ie)yp^ 

we see from the proof of Proposition 14 that under our assumptions 

(282) lim E(Fj^^(e*^)) = 

N^oo ’ ’ 

In fact, this convergence is uniform in 6 as the relevant continuum Toeplitz 
determinant is independent of 9 due to the fact that the law of the eigen¬ 
values of the CUE is rotation invariant. 
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Due to the convergence in distribution of Proposition 14, there exists a 
probability space where one can construct a sequence of random variables 

and (-^i, -^l) such that for each N,M 

/ N N \ 

(283) 

Vfc=i k=i J 

where {zi,...,zn) is sampled from P 7 v,M; are independent complex 

Ganssians with real and imaginary parts independent centered real Gaus- 
sians of variance j'/2, and 

(284) ^ 
almost surely as —)• oo and M — N —)• oo. 

This implies that if we construct an object agreeing in law with ^ 
from these qnantities, it converges uniformly almost surely to 


(285) 


^-/3ReEt=i 


Due to this uniform convergence (as well as the fact that {—Zj)j = {Zj)j) 
and the convergence of E(F^^^(e*^)) to it follows that the cor¬ 

responding object whose law agrees with that of ^ ^ converges almost 

surely to as —)• oo in such a way that M — N ^ oo (the convergence is 
in the topology of weak convergence of measures). Thus we have the claim. 

□ 


6.2. Proof of Proposition 18, Proposition 17, and Proposition 16. 

We begin with Proposition 18 

Proof of Proposition 18. This follows directly from combining Theorem 26, 
Proposition 30, and Lemma 49. □ 

We then move onto Proposition 17. 

Proof of Proposition 17. This is essentially just a combination of onr esti¬ 
mates on the Fredholm determinant and the strong Szego theorem or results 
on the asymptotics of Toeplitz determinants with Fisher-Hartwig singular¬ 
ities - namely we combine Theorem 25, Theorem 26, Proposition 30, and 
Lemmas 48, 49, and 50. As we make statements about uniformity, we give 
a few more details now. 

1) Gonvergence follows directly from Proposition 30, Theorem 25, and 
Lemma 48. Uniformity follows from the fact that the continunm Toeplitz 
determinant is independent of 6 (e.g. by the Heine-Szego identity and the 
fact that the law of the eigenvalues of a GUE matrix is rotation invariant), 
while onr estimate in Lemma 48 is nniform. 

2) This is as the previous case except we use statement 1) of Theorem 
26 instead of Theorem 25 for the asymptotics of the continuum Toeplitz 
determinant. 
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3) The reasoning is the same again - Proposition 30 gives the factoriza¬ 

tion into a Fredholm determinant with uniform estimates and a continuum 
Toeplitz determinant, which we’ve denoted by E^,Lid,S'). The fact that 
this increases to is proven for example in [31]. 

4) Again we make use of Proposition 30, Theorem 26, and Lemma 49. 
The uniformity now follows from the fact that the continuum Toeplitz de¬ 
terminant is rotation invariant (so we can rotate 9 from the non-singular 
part to 6' — 6 in the singular part) and make use of the uniformity of state¬ 
ment 1) of Theorem 26 (now there is only one singularity so the estimate is 
uniform in its location). 

5) and 6) Follow directly from combining Proposition 30, Theorem 26, 
and Lemma 49 and Lemma 50 respectively. 

□ 


Finally we prove our variance estimate. 

Proof of Proposition 16. Let us write down explicitly all the definitions. We 
first expand the square and then note in each term we can interchange the 
order of integration as everything is non-negative. 


E 


f27r 




= r r I,MM' 

Jo Jo "" ® E(FL.r(e*^))E(Ff,,,,(e*®')) 27 r 27 r 

r27r /* 27 r 

(286) - 2 / / /(e*®)/(e*®') 

Jo Jo 

r27T p27t 

+ / f(e'^)f(E^') _ 

Jo Jo E(Fj,^(e*e))E(F^,^(e*e')) ^vr 




dO dO' 




dO dO' 


=■ h — 2/2 + h- 


Combining parts 1) and 3) of Proposition 17, we see that (even if we just 
assume M — N^ooasN^ 00 ), 


p27z i'2'k 

hm h= / fie'^)fiE^') 
Jq Jq 

(287) = [ [ 

Jo Jo 


E(e/3W(0)e/3W(e')) 


]E(g/3Xi,(0))]E('g/3Xi(0')) 2 tt 2 -k 

0 Jo 27r 27r 

Combining parts 1), 2), and 4) of Proposition 17, we have (as N/M 


0) 


TV^oo " io io 27r27r 
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For I 3 , we repeat an argument from [6]. We split the domain of integration 
into \9 — 9'\ > 2to and \9 — 9'\ < 2tQ. For the first case, we use 2) and 5) of 
Proposition 17 and find 


lim / 

J\e-e'\>2to 




w 




d9 d9' 


(288) 


'|6»-6»'|>2io 


ZTT ZTT 


For the latter case, we note that if N\9 — 9\ < 1, we note that by the 
asymptotics of a and the fact that log(a;“^ sin(a:/2)) is bounded for small x. 


(289) E(F^,,(e*^)F^^^(e*^')) = 0(l)fV^^ 

where the 0(1) is uniform in A^|0 — 0'| < 1, so we see that 


(290) 


f E«M(^")4M(e")) 

ie-9'\<N-i ^E(F^^^(e*'^))E(F^^^(e*®')) 27r 27r 

= 0(7V-1)0(1)^ 

N— 

= 0 {N^-^) 


which tends to zero as —)■ 00 under our assumptions. 

On the other hand, for the N~^ < \9 — 9'\ < 2to case we write 


E (<M(e")<M(e"')) 


0 (l)e^ 





(291) 


X e 


j--iN\e-e'\ 



log|0-0'| 


where the 0(1) can be taken to depend only on /? (so there are precise 
asymptotics for it, but we care only about it being uniformly bounded for 
hxed /3). Here we made use of the fact that |2sin(x/2)| = |e** — 1|. Again 
making use of the s —)• 0 asymptotics of a the hrst integral is hnite (and its 
value only depends on /3), while the second integral we write as 


(292) 


.-iN\e-9'\ ^ 02 o2 

- ds —— log N —— log \9 — 9'\. 


Making use of the s —>■ 00 asymptotics of u, we see that the above integral 
is uniformly bounded and we conclude that 
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^E«^(e*^))E(4^(e*^'))2vr2vr 

“ y|e-e'|<2io 27r27r 

As limjv^oo -^3 must be independent of to, and this tends to zero as to 0 
(since /3^ < 2), putting everything together we see that 


f*27r r27T 

(294) hm h = j / /(e*®)/(e*®')| 

Jo Jo 

We conclude that for N/M —>■ oo as —)• oo, 




27r 27r 


lim E 

N^oo 


r' 27 r 




(295) 


/* 27 r /* 27 r 

/ / /(e*^)/(e*^') 

Jo Jo 




27r 27r 


As 


(296) lim 

L^oo 

Fatou’s Lemma implies that 


Jd 



je 


^ZTT pZlT 

/ / /(e‘'')/(e“‘')|e-- e-I » 

Jo Jo 


*61'I-4^ dO dO' 
2 tt 27r 


(297) <limsup f " f ^ 

L^oo Jo Jo 27r 27r 

but (295) implies the opposite relation (the left side of (295) is non-negative 
- it’s a variance) so we see that 


(298) 


lim lim E 
L—>-oo N—>-oo 



^)if^N,M,Lidd) - A,MidO))) 


= 0 . 


□ 


6.3. Proof of Theorem 9. We now prove our main result. 

Proof of Theorem 9. We wish to prove that for any continuous / : T —)■ 
[0, oo), as A^ —>• oo (and N/M —>• 0) 

r2'K r27r 

(299) / f{e^VNMide)A f{PV{de). 

Jo Jo 

We write 



A DISCRETE LOG GAS, DISCRETE FH-SINGULARITIES, AND GMC 


59 


/*27r /*27r 

/ f{e^VNM(de)= 

r2n 

(300) +/ f{e^^){AMidO)-AMLidO)). 

Jo 

By Proposition 13 and the definition of /i^, 

r27r 

(301) / f{e^VN,M,LidO) ^ / f{eVm 

Jo Jo 

if we first let —>• oo and then L —> oo. 

On the other hand, as we saw in Proposition 16 that 


(^(/ /(^*^)(/^^,m(^'^)-/^w,m,l(^^))) ^ = 0 , 

we have in particular that if we first let —)• oo and then L —)• oo, 

(303) / f{e^^){fi^]^{d6) — ^ 0. 

J 0 

We conclude, by Slutsky’s theorem, that if NjM —)• 0 as A^ —)• oo, then 

/•27r /•27r 

(304) / /(e*^)/,J^(d0)4 / fie^VidO) 

Jo Jo 

as A^ —)• oo. 

□ 


7. Discussion 

We now discuss less precisely what might happen to the Fredholm deter¬ 
minant det(/ -|- K) in the case where N/M —>■ (7 S (0,1) as A^ —^ 00 . For 
simplicity, we focus on the case where the distance between the singulari¬ 
ties is bounded away from zero. From our analysis, it looks reasonable that 
tr(it''^) = 0{{N/My) so if we assume that q is small, to argue that the 
asymptotics of det(/ + K) differ from one, it should be enough to argue that 
tr(it') 0. 

Looking at our argument, we see that the most relevant contribution to 
the trace comes from values of the kernel between [M“^,Ai“^] (close to 
the point M~^ - elsewhere there seems to be a term of the form e“^/'^). 
Working a bit harder on the asymptotics of the orthogonal polynomials, one 
could check whether or not there is some cancellation in this regime (e.g. 
is there a symmetry causing the integrals along (1 — A^“^, 1 — M~^) x Wj 
and (1 -|- M~^, 1 -|- A^“^) x Wj to cancel). As the analysis of the asymptotics 
was a bit heavy as it is, we choose not to pursue this question further here, 
but note that it seems unlikely that there would be such a cancellation 
to all orders in q, which implies that it is quite reasonable to expect that 
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the Fisher-Hartwig formula is unlikely to hold as it is for q > 0. It would 
be an interesting problem to try to calculate more precise asymptotics of 
the Fredholm determinant in the N/M —)• 5 > 0 case. Another interesting 
questions is how this affects the law of the limiting chaos measure or perhaps 
the law of the maximum of the characteristic polynomial (for more on this 
type of questions in the context of characteristic polynomials of log-gases, 
see [14, 15]) 

We wish to underline the interesting contrast between this and the be¬ 
havior of the linear statistics. If one just has M — N^ooasN^ 00 , the 
law of any finite collection of the linear statistics agrees with the continuum 
case, so the linear statistics are essentially unable to see the discreteness of 
the model, while the characteristic polynomial does seem to see it. This 
is presumably due to the fact that the characteristic polynomial essentially 
contains information of all polynomial linear statistics, so it might be reason¬ 
able to expect that the polynomial linear statistics also see the discreteness 
when one takes the power to grow with N and M. 
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